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In these lectures we discuss some elementary concepts in connection with the theory of 
symmetric spaces applied to ensembles of random matrices. We review how the relation- 
's! \ ship between random matrix theory and symmetric spaces can be used in some physical 
^ , contexts. 
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1 Introduction 



1.1 Outline of the lectures 

In these lectures, that have been elaborated in part from the review by the author and 
Caselle , we will deal with some elementary concepts relating to the description of random 
matrix ensembles as symmetric coset spaces. Some figures have been added to the text 
as illustrations. These have been borrowed from various reviews and papers on random 
matrix theory with permission from the authors and publishers. 

As the conference organizers asked me to give elementary lectures, these lecture notes 
should be accessible to a wide audience. Even though none of the material presented here 
is entirely new, perhaps these lectures can still serve the purpose of stimulating interest in 
the research on random matrices. 

In the first lecture we will introduce some of the physical systems where Random Matrix 
Theory (RMT) is a useful tool. Then we will give a brief overview of the most important 
elements in random matrix theory. Part of the material here is based on excerpts from 
the excellent review by Guhr, Muller-Groeling and Weidenmiiller [2j. Since the audience 
consists of mathematicians, most of whom are not working in this field, we will try to be 
as clear as possible in describing what RMT is used for, why it is relevant, and how its 
predictions are compared to experimentally or numerically measured spectral fluctuations. 
In particular we will discuss the concepts of unfolding and universality. 

The second and third lectures are essentially a shortened version of some of the material 
presented in P, with some modifications and additions. In the second lecture we will see 
that hermitean random matrix ensembles can be identified with symmetric coset spaces 
related to a compact symmetric subgroup. This was realized early on by Dyson (3J 0] and 
by Huffmann j^j, but the advantages due to this fact for a long time was not frequently 
understood by physicists in applications to the above-mentioned systems. The well-known 
theory of simple Lie algebras and symmetric spaces [HI IE] can be applied in random 
matrix theory to obtain new results in various physical contexts pp. Key elements are 
the classification of symmetric spaces in terms of root systems and the theory of invariant 
operators on the symmetric manifolds. 

We will discuss elements from the theory of symmetric spaces, using explicit examples 
drawn from low-dimensional Lie algebras. As we will see, symmetric spaces of positive, 
zero, and negative curvature correspond to well-defined types of random matrix ensembles. 
Along the way, we will also discuss the identification of various elements of RMT with 
the corresponding quantity on the symmetric space. In particular it will be clear that 
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eigenvalue correlations in random matrix theory has a geometric origin in the root systems 
characterizing the symmetric space manifolds. 

In the third and last lecture we discuss a few examples of applications of symmetric coset 
spaces in random matrix theory. A few more topics have been discussed in £Q . As the first 
and most important application, we discuss the classification of disordered systems arising 
from the Cartan classification of symmetric coset spaces. Most of the hermitean random 
matrix ensembles corresponding to symmetric spaces have known physical applications. 
We will not have time to discuss these here, but physical applications of the ensembles 
appearing in Table El along with more random matrix ensembles, were discussed or at 
least mentioned in £Q , where references to the literature can also be found. (I apologize to 
those authors whose work has been overlooked here, as I am certain not to be aware of all 
the published work on applications of random matrix ensembles.) 

Our second example is the solution of the DMPK equation using the theory of zonal 
spherical functions (the eigenfunctions of the radial part of the Laplace-Beltrami operator 
on the symmetric space). The DMPK equation is the scaling equation determining the 
probability distribution of the transmission eigenvalues for a quantum wire as a function 
of the length of the wire. It will be identified with the equation of free diffusion on the 
symmetric space. 

The last example we give is more weakly related to RMT (it is related only through 
the diffusion equation) and concerns the integrability of certain Id models referred to as 
Calogero-Sutherland models. There is a connection between these models and the theory 
of Lie algebras: Calogero-Sutherland models are closely related to root systems of Lie 
algebras or symmetric spaces. Olshanetsky and Perelomov [0] provide an exact statement 
as to when these models are integrable and directly express the physical integrals of motion 
as the algebra of Laplace operators related to the Lie algebra or symmetric space. These 
results lead to a detailed list of spectra and wave functions for a variety of quantum systems. 



1.2 A few comments on the classification scheme 

It is interesting to note that a wide variety of microscopically different physical systems can 
be described by the same type of spectral fluctuations. RMT describes the generic features 
of the spectrum, without regard to dynamical principles or details of the interactions. This 
allows a separation of generic spectral fluctuations from system-specific ones. The only 
input in RMT is symmetry through the postulated invariance of the various ensembles. 
This scenario leads to a classification of physical disordered systems into symmetry classes 
characterized by universal spectral behavior. 
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In the context of the classification scheme a new paper by Heinzner, Huckleberry and Zirn- 
bauer JU] should be mentioned, even though none of the material there will be covered 
in these lectures. In this paper, the authors prove the correspondence between symmetric 
spaces and symmetry classes of disordered fermion systems with quadratic Hamiltonians. 
This is done by considering both unitary and antiunitary symmetries and then removing 
the unitary symmetries by considering the decomposition of the space of good Hamiltoni- 
ans into blocks associated with unitary subrepresentations in the Nambu space of fermionic 
field operators. The relevant structures on this space are transferred to a space of homo- 
morphisms where the unitary symmetry group acts trivially. The various cases occurring 
for the remaining anti-unitary symmetries then lead to the classification in the symmetric 
space picture. The authors also observe that when second quantization is undone (as it is 
in the physical systems corresponding to the new symmetry classes arising in physics in 
addition to Dyson's symmetry classes), a remnant of the canonical anticommutation rela- 
tions of the fermionic field operators is imposed on the Nambu space. This is the reason 
why we get new structures in the physical context of disordered fermions. 



1.3 Topics outside the scope of these lectures 

Of course, a significant fraction of the (vast) literature on the theory and applications of 
RMT will not be covered here. This applies also to various types of extensions of simple 
hermitean random matrix models, multimatrix models, and several phenomena in random 
matrix theory, e.g. parametric correlations and phase transitions, and theoretical issues 
like universality proofs and the supersymmetric formalism. The obvious reason is that our 
focus should be on the relationship with symmetric coset manifolds. A more extensive 
introduction to random matrix theory as well as experiments and numerical simulations 
performed to disclose RMT behavior in spectra can be found in [2], that gives a good 
overview of the various aspects of random matrix theory. 

Note that even the list of simple hermitean random matrix ensembles mentioned in these 
lectures is not complete, though we do discuss the main ones. A few more ensembles were 
briefly discussed in p. In addition there is a large number of non-hermitean random matrix 
ensembles. There has been some recent activity in this field, where most of the papers have 
dealt with the problem of finding the eigenvalue distribution in the complex plane. The 
concept of orthogonal polynomials has also been extended to the non-hermitean case (see 
e.g. HH). 

In the applications discussed here, random matrices are used to describe statistical fluc- 
tuations in the spectra of quantum operators. In some field theoretical contexts, random 
matrices are employed in a different way. In quantum gravity, they describe random dis- 
cretized surfaces corresponding to string world sheets of different genus. The partition 
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function corresponding to an integral over the gravitational field is the sum of such sur- 
faces, much like in the large N expansion in QCD due to 't Hooft. In such a context, it is 
the field rather than a quantum operator that is substituted by a random matrix. We will 
not discuss these applications of random matrices. For an introduction see [T2] . 



2 Chaotic systems and random matrices 

The wide range of physical problems where random matrix theory can be successfully 
applied has made it into an important branch of mathematical physics. As an introduction 
to the topic, we will give a brief historical survey of the developments which led to the 
present situation. 



2.1 Many— body systems and Wigner— Dyson ensembles 

In the 1950's, Wigner ^H] developed a theory of random matrices to deal with resonance 
spectra of heavy nuclei. Experiments with neutron and proton scattering gave precise 
information on levels far above the ground state, whereas nuclear structure models could 
only predict the positions of levels close to the ground state. Wigner conceived of a 
new way of studying the spectrum: a statistical theory that could not predict individual 
energy levels, but described, in Dyson's words, "the general appearance and the degree of 
irregularity of the level structure" [Sj. This provided a tool for studying complex spectra. 
Wigner's theory dealt with ensembles of random matrices modelling the Hamiltonians of 
nuclei. In the early 1960's, Dyson developed Wigner's approach further in a series of 
papers [3] where he treated scattering matrix ensembles. Typical of the spectra obeying 
Wigner-Dyson statistics is that energy levels are correlated and repel each other. Such a 
level repulsion is not present in spectra obeying Poisson statistics (cf. Fig. QJ. 



Shortly thereafter RMT was applied to the spectra of atoms. Later on, modern laser 
spectroscopy has allowed a comparison with the complex spectra of polyatomic molecules. 
Nuclei, atoms and molecules are all examples of complex many-body systems with a large 
number of degrees of freedom. 
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Figure 1: Nearest neighbor spacing distribution for the "Nuclear Data Ensemble" com- 
prising 1726 spacings (histogram) versus s = S/D with D the mean level spacing and 
S the actual spacing. Note the level repulsion in random matrix theory (GOE denotes 
the Gaussian Orthogonal random matrix Ensemble) at zero spacing as compared to the 
Poisson distribution. Taken from reference Used with permission. 
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Figure 2: Four important billiards: the Bunimovich stadium (top left), the Sinai billiard 
(top right), the Pascalian snail (bottom left) and the annular billiard (bottom right). 
Reprinted from reference [2] with permission from authors and from Elsevier. 

2.2 Quantum chaos 

It was later realized that the Wigner-Dyson ensembles could be applied to the description 
of chaotic quantum systems. In these, a quantum particle is reflected elastically at the 
boundaries of some given domain, so that the total energy is constant. In such a system, 
all constants of the motion (except the energy) are destroyed by randomness, and there 
are no stable periodic orbits. The system is referred to as a quantum billiard (see Fig- EJ) - 
Due to the shape of the domain, normally chosen to be two-dimensional, the trajectory of 
the particle is completely random. Such a system may have just a few degrees of freedom 
and may be simulated by microwave cavities. The reason is that the wave equation for 
the electromagnetic field in the cavity has the same form as the Schrodinger equation 
for a two-dimensional quantum billiard, if the geometry of the cavity and the boundary 
conditions are chosen appropriately. In this Schrodinger equation, the Hamiltonian for the 
free particle is simply the Laplace operator. With appropriate boundary conditions the 
system may also be equivalent to a vibrating membrane. 

Spectra of such systems with up to a thousand eigenmodes have been studied experimen- 
tally (for a list of references, see [2j). The results can be interpreted in terms of classical 
chaos. The general picture emerging from such experiments is that if the shape of the 
cavity is such that the classical motion in it is integrable (regular), the spectrum behaves 
according to a Poisson distribution. If the corresponding classical motion is chaotic, the 
spectral fluctuations behave in accordance with the Wigner-Dyson ensembles. This is also 
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the content of the famous (but unproven) Bohigas conjecture. 



2.3 Mesoscopic systems, BdG and transfer matrix ensembles 

A newer class of interesting quantum systems with randomness is provided by disordered 
metals, i.e. micrometer-sized metal grains, of which we wish to study the transport prop- 
erties when the grain is connected to electron reservoirs through (ideal) leads (for a review, 
see [X5J ) - Such microstructures can be fabricated in a cavity in a semiconductor. They 
can be made so thin that the electron gas inside them is effectively two-dimensional. Such 
structures are quantum systems, yet they are large enough for a statistical description to 
be meaningful. Therefore they are referred to as mesoscopic systems. The transmission 
eigenvalues related to the transfer matrix formalism determine the conductance, which is 
the main observable. We will discuss these systems more explicitly in the third lecture. 

Disorder in such systems arises because of randomly distributed impurities in the metal. 
The electrons in the conduction band are scattered elastically by the resulting random 
potential as we apply a voltage drop across the metal sample. Such a system is called 
diffusive. At low temperature (below IK), inelastic electron-phonon scattering, which 
changes the phases of the electrons in a random way, can be neglected. The electrons are 
therefore phase coherent over the length of the sample. In the diffusive regime (with the 
conductance decreasing linearly with sample size) we expect Wigner-Dyson statistics to 
apply to the spectrum of energy levels (at least up to an energy separation given by the 
so-called Thouless energy). The random matrix theory of quantum transport, however, 
deals mostly with the transmission eigenvalues. It relates the universality of transport 
properties (i.e., the independence of these of sample size, degree of disorder, etc.) to the 
universality (to be discussed) of correlation functions in random matrix theory. 

The behavior of the conductance depends critically on the dimensionality of the sample. 
In Id all states are localized and the sample is an insulator, while d = 2 is the critical 
dimension and derealization occurs for d > 2. If the sample has the shape of a grain, it is 
called a quantum dot. A quantum wire is a very thin (quasi-one dimensional) metal wire 
with similar properties, that allows to study the scaling properties of observables related 
to transport as a function of the length of the wire through a generalized diffusion equa- 
tion. Quasi-ld wires are particularly good laboratories for testing RMT. Since they are 
not strictly one-dimensional, they possess a diffusive regime (cf. the comment above on di- 
mensionality); still, non-perturbative analytical methods are applicable (scaling equation, 
field theoretical methods). 

The phenomenon of localization occurs when the previously extended Bloch wave functions 
of the multiply scattered electrons are cancelled by destructive interference. This happens 
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Figure 3: Electron micrograph of a quantum dot in the shape of a Bunimovich stadium, 
with 1 [xm bar for scale. The electrons can move in the black area. Two leads are coupled 
to the stadium. Reprinted with permission from ref. [T^j. Copyright 2005 by the American 
Physical Society. 

when the density of impurities reaches a critical value. As a consequence, the metal sample 
goes from being a conductor to being an insulator. In the localized (insulating) regime 
the conductance decays exponentially as a function of sample size and the typical length 
scale of the decay of electron wave functions is given by the localization length £. In the 
localized regime, wave functions of nearly degenerate states may have a very small overlap. 
As a consequence, the level repulsion typical of the Wigner-Dyson ensembles is suppressed 
and we expect Poisson statistics (uncorrelated eigenvalues) on length scales larger than £. 

We may also consider a metal grain in which electron scattering takes place at the bound- 
aries. In the latter case the mean free path of the electrons is large compared to the size 
of the system and we speak of a ballistic system. A ballistic quantum dot is very similar 
to a billiard and is described by Wigner-Dyson statistics (see Fig. |3J). 

In addition to these systems, heterostructures consisting of superconductors in conjunction 
with normal metals are successful candidates for a random matrix theory description. In 
these we have particle-hole symmetry and the (first-quantized) Hamiltonians in this case 
are described by four new so-called Bogoliubov-de Gennes (BdG) ensembles [T7] . 

The existence of several mesoscopic regimes (localized, diffusive, ballistic, with varying 
degree of disorder) enriches the phenomenology and extends the applications of RMT 
beyond the ones discussed for many-body and chaotic systems. In particular, the case 
of normal-superconducting quantum dots leads to four entirely new symmetry classes of 
Hamiltonians. 
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In mesoscopic systems, the operator that is modelled by random matrix ensembles is either 
the Hamiltonian or the transfer matrix. In the random Hamiltonian approach, ensemble 
averages are calculated in the supersymmetric formalism ^H] using scattering theory. (The 
super symmetric formalism is partly discussed also in j2J - In addition see fH] for QCD- 
related applications.) In the random transfer matrix approach, which we will discuss more 
explicitly later, the transmission eigenvalues determine important transport properties, of 
which the most central one is the conductance determined in the Landauer-Lee-Fisher 
theory. To every ensemble of Hamiltonians corresponds an ensemble of transfer matrices 
(and one of scattering matrices as well) describing the same physical system. 



2.4 Field theory and chiral ensembles 

Random matrices are also useful in relativistic quantum field theory. In a gauge field 
theory, the interactions between the gauge field and the fermions is expressed as a gauge 
field integration over the fermion determinant, obtained by integration over the fermionic 
(Grassmann) variables in the partition function. The fermion determinant involves the 
Dirac operator and depends on the gauge fields. This operator can also be modelled 
by an ensemble of random matrices. If chiral symmetry is present, these ensembles are 
called chiral random matrix ensembles. The matrices in these ensembles have a block- 
structure similar to the Dirac operator in the chiral basis, and may possess zero modes. 
This application of random matrix theory was developed by Verbaarschot et al. in the 
1990's (for some early work, see |20j). 

The procedure of substituting the Dirac operator with a random matrix makes it possible to 
perform the integration over the gauge fields by integrating over the ensemble of random 
matrices. By using techniques borrowed from the supersymmetric formalism originally 
developed for scattering problems in the physics of quantum transport, one can then obtain 
the effective low-energy partition function in the gauge theory [25- In the presence of 
spontaneous symmetry breaking, it is expressed as an integral over the Goldstone manifold 
(i.e. over the degrees of freedom associated with spontaneous symmetry breaking in the 
gauge theory). This way the symmetry breaking pattern is also obtained and in addition, 
sum rules for the eigenvalues of the Dirac operator can be derived. Recently, there have 
been interesting new developments in this field, relating the QCD partition function to 
integrable hierachies (for a comprehensive review with references, see the lecture series by 
Verbaarschot |2*2*]). 

The order parameter for the spontaneously broken phase (the quark condensate in QCD) 
is proportional to the density of Dirac eigenvalues at the origin of the spectrum (i.e., at 
zero eigenvalue). This is expressed by the Banks-Casher formula. Therefore the Dirac 
spectrum, and especially the regime around the origin, is of major interest. Its generic 
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features can be studied in RMT. A wealth of numerical simulations [231 * n lattice gauge 
theory show that data agree with RMT predictions for level correlators and low-energy 
sum rules for the Dirac eigenvalues. This shows that fluctuations of Dirac eigenvalues 
are generic and independent of details of the gauge field interactions. Two reviews with 
references to original work are given in . 

Note that since these lectures are introductory, we will limit the discussion to hermitean 
random matrix models. Therefore we will not discuss the new developments in this field 
involving non-hermitean random matrices (see e.g. |2S|), even though these are very 
intriguing as they allow a study of the case of nonzero baryon chemical potential. The 
same is true for the new research developments in non-hermitean random matrix theory 
in other branches of physics. Such applications include for instance the description of 
non-equilibrium processes, quantum chaotic systems with an imaginary vector potential 
(related to superconductivity), and neural networks. 



3 What is random matrix theory? 



We have established that some of the things that may make a system chaotic are: complex 
many-body or gauge field interactions, a random impurity potential, or irregular shape 
(where irregular means any shape leading to classical chaotic motion). Because of the high 
number of degrees of freedom or the complexity of the interactions in chaotic systems, the 
quantum mechanical operators whose spectra and eigenfunctions we are interested in may 
be unknown or inaccessible to direct calculation. Such an operator may be a Hamiltonian, 
a scattering or transfer matrix, or a Dirac operator, as we have seen above. These operators 
can be represented quantum mechanically by matrices acting on a Hilbert space of states. 

In studying nuclear resonances, Wigner ^2! proposed replacing the Hamiltonian of the 
nucleus by a random matrix taken from some well-defined ensemble. How do we determine 
which ensemble is appropriate? Wigner proposed choosing it so that the general global 
symmetries of the underlying physical problem are preserved. This means that we should 
choose random matrices with the same global symmetries as the physical operator we are 
modelling. In doing so, we will preserve the spectral properties that depend only on these 
symmetries. These are the aspects of the spectrum we will be able to study using RMT. 

As it turns out, the gross features of the random matrix theory spectrum, like for instance 
the macroscopic shape of the eigenvalue density, is usually not of interest for studying 
the physical system, as it does not even remotely resemble the actual spectrum and is 
non-universal in that it depends on the choice of random matrix potential. The strength 
of the RMT description lies in the universal behavior of the eigenvalue correlators in a 



12 



certain scaling limit. As we will see, this leads to universal spectral fluctuations that are 
reproduced by a multitude of physical systems in the real world, several of which we have 
discussed in the previous paragraph. 

The statistical theory of Wigner ^3] and Dyson jH] is a theory of random matrices. These 
are matrices with random elements chosen from some given distribution. In the applications 
to be discussed here, we use them to represent a quantum mechanical operator, whose 
eigenvalue spectrum and eigenstates are of physical interest. In the classical example of 
the heavy nucleus, the Schrodinger equation 

HV i = E i % (3.1) 

gives the energy spectrum of the nucleus with Hamiltonian H. In case this Hamiltonian 
is unknown or too complicated, we may choose to study the eigenvalues of an appropriate 
ensemble of random matrices instead. The ensembles studied by Dyson in the early sixties 
were ensembles of scattering matrices, but the principle is the same. 

The random matrix eigenvalues should behave much like the energy levels of the nucleus, 
if we choose our ensemble of random matrices appropriately. To achieve this we analyze 
the physical symmetries of the Hamiltonian and choose an ensemble of random matrices 
with these same symmetries. The characteristics of the spectrum that depend only on 
these symmetries will be present also in the spectrum of the random matrix eigenvalues. 
These are the universal features of the system, and they are obtained after removal of 
the non-universal behavior through a proper rescaling procedure. They do not depend on 
details of the physical interactions or, in the RMT, on the choice of RMT potential. 

The number of fields where random matrix models are applied nowadays is growing fast. 
To convey the main ideas, we will concentrate on a few major applications. 



4 Choosing an ensemble 

What are the physical symmetries that determine the ensemble of random matrices? 
Dyson's analysis shows that for the Hamiltonian ensembles there are three generic cases: 

1) The system is invariant under time-reversal (TR) and the total spin is integer (even- 
spin case), or the system is invariant under time-reversal and space-rotation (SR) with no 
regard to the spin. 1 

1 Time reversal invariance arises from the fact that (in the absence of magnetic fields), if ^(r,t) is a 
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2) The system is invariant under time-reversal but with no rotational symmetry and the 
total spin is half-odd integer (odd-spin case). 

3) The system is not time-reversal invariant. 

In all these cases, the ensemble of random matrices is invariant under the automorphism 

H -> UHU' 1 (4.1) 

where U is an orthogonal (case 1), unitary (case 3) or symplectic (case 2) matrix. It is 
this property that gives rise to the identification of the manifolds of random matrices with 
symmetric spaces. This identification is very useful, since it leads to the possibility of 
applying the known theory of symmetric spaces in physical contexts. 

The properties of these ensembles are summarized in Tabled 2 

Table 1: Dyson's three universality classes of random matrices. The Dyson index (3 counts 
the number of degrees of freedom of the matrix elements (each matrix element is a real, 
complex or quaternion number) and is used to specify the ensemble. The ensembles are 
named the orthogonal, unitary or symplectic ensemble, respectively, according to the sub- 
group of invariance. 



p 


TR 


SR 


H 


U 


1 


X 


X 


real, orthogonal 


orthogonal 


2 




(X) 


complex, hermitean 


unitary 


4 


X 




real quaternion, self-dual 


symplectic 



The relevant physical symmetries depend on the system under consideration. Novel en- 
sembles not included in this classification arise if we impose additional symmetries or 

solution of the Schrodinger equation, so is \&*(r, — t). We define the action of the (anti-unitary) time- 
reversal operator T on a state $ by T^S = K^f*, where if is a unitary operator. T should reverse the 
sign of spin and total angular momentum. This requirement can be satisfied by the choice K = e" s y for 
spin rotation or K = e l7Tjy for space rotation, with a standard representation of spin or space rotation 
matrices. Now T 2 = KK* = e 2lnSy = (— )" = ±1. These two cases correspond to integer and half-odd 
integer spin (or presence or absence of space rotation invariance in case of J y ). In the former case, K can 
be brought to unity by a unitary transformation $ — ► U^S of the states, during which K transforms into 
UKU T . Once this choice has been made, the only transformations on H and $ allowed are H — > RHBT 1 , 
$ — > Rty with R an orthogonal matrix. In the latter case, a block-diagonal form for K may be chosen 
and only symplectic transformations are allowed on H and ^ . In case there is no time reversal invariance, 
unitary transformations on H and are allowed. 

2 The dual Q R of a matrix Q consisting of N x TV real quaternions q = qo + q • r is Q R = q~ji where 
q = qo — q • r and n = —ioi where Oi are Pauli matrices. Self-dual means Q = Q R . 
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constraints. This is the case for the four universality classes of the Hamiltonians of normal- 
superconducting (NS) quantum dots (for a detailed discussion see [HI)- Also here the four 
universality classes are distinguished by the presence or absence of TR and SR, but with 
the difference that the Hamiltonian possesses so-called particle-hole symmetry. As elec- 
trons tunnel through the potential barrier at the NS interface, a Cooper pair is added to, 
or subtracted from, the superconducting condensate. Equivalently, an electron incident on 
the interface is reflected as a hole (a phenomenon referred to as Andreev reflection). 

Three more chiral symmetry classes are realized in chiral gauge theories. In this case the 
ensemble is chosen such that the chiral symmetry [75, %lp\ = 0, the zero modes and possible 
anti-unitary symmetries [Q, % p\ — of the Euclidean Dirac operator (where Q is some 
anti-unitary operator) are reproduced by the ensemble. The properties of Q determine if 
there is a basis in which the Dirac operator is real or quaternion real. In this case we also 
have three symmetry classes distinguished by the index j3. 

As we will see later, to every hamiltonian ensemble there corresponds a scattering matrix 
ensemble and a transfer matrix ensemble. The hamiltonian and scattering matrix ensem- 
bles of a given physical system are related to each other in a simple way: the hamiltonian 
ensemble is simply the algebra subspace corresponding to the compact symmetric space 
of the scattering matrix. Because of the way transfer matrices are defined, the ensemble 
of the transfer matrix of the same system is not the corresponding non-compact space. 
Nevertheless, it will neatly fit into the classification scheme. For scattering and transfer 
matrices in the theory of quantum transport we use ensembles constrained by the physical 
requirements of flux conservation, time-reversal symmetry, and spin-rotation. Flux con- 
servation imposes the condition of unitarity on the scattering matrix, and determines the 
symmetry class of the transfer matrix too. We will discuss some explicit examples in the 
third lecture. 



5 General definition of a matrix model 



Since many of the interesting quantum operators are physical observables, we often deal 
with hermitean random matrices. Non-hermitean operators are also of interest. The 
number of fields where they are important is growing, and so is the research effort devoted to 
this large class of random matrix theories 3 . Here we will only be concerned with hermitean 

3 Non-hermitean quantum operators are applied in non-equilibrium processes and dissipative systems. 
Examples are quantum chaotic scattering in open systems, conductors with non-hcrmitcan quantum me- 
chanics, and systems in classical statistical mechanics with non-hermitean Hamiltonians. They are also 
useful in schematic random matrix models of the QCD vacuum at non-zero temperature and/or large 
baryon number, where the fermion determinant becomes complex. 
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random matrices. 



A hermitean matrix model is defined by a partition function 



Z ~ J dHP(H) (5.1) 

where if is a square N x N hermitean matrix with random elements Hij chosen from 
some given distribution. At the end of our calculation, we will take the limit iV — > oo to 
counteract the fact that we are using, for technical reasons, a truncated Hilbert space. In 
equation ()5.1|) P(H) is a probability distribution in the space of random matrices and dH 
is an invariant measure in this space. Such a measure is required to give physical meaning 
to the concept of probability: P(H)dH is the probability that a quantum operator in 
the ensemble will belong to the volume element dH in the neighborhood of H. Since the 
ensemble of matrices {H} does not in general form a group, the existence of such a measure 
is not entirely trivial. Dyson [3| showed that a unique invariant (Haar) measure dH exists 
for the three ensembles discussed in the previous section. 

One can show that a probability distribution of the form 

P(H) ~ e - ctrV{H) (5.2) 

is appropriate (see j2E] for more details on this point). Note that such a weight function 
is needed to keep the integrals from diverging. Here c is a constant and V(H) is a matrix 
potential, typically a polynomial with a finite number of terms. The simplest choice is a 
quadratic potential, in which case the matrix model is called gaussian. A common choice 
for the Wigner-Dyson ensembles is 

Pp{H) ~ e"f> H2 (5.3) 

where v has the same dimension as the matrix elements of H. The factor in front of the 
potential is chosen proportional to j3 for later convenience and proportional to N so that the 
spectrum, that has support on some finite interval on the real axis, will remain bounded in 
the large N limit. In this case the macroscopic eigenvalue spectrum is a semicircle. Such a 
spectrum is quite unrealistic for most physical systems, and not interesting in itself. What's 
more, the macroscopic form of the spectrum depends on the form of the potential. However, 
the statistical eigenvalue fluctuations will turn out to be universal and independent of the 
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matrix potential in the large- N limit, if we scale the eigenvalues appropriately (a process 
referred to as unfolding) . Our focus will be on such universal quantities. 



The probability P(H)dH is invariant under the automorphism 



H -> U- l HU (5.4) 

of the ensemble to itself, where U is a unitary N x N matrix. By doing an appropriate 
similarity transformation (|5.4jl on the ensemble of random matrices, the Haar measure dH 
and potential V(H) can be expressed in terms of the eigenvectors and eigenvalues of the 
matrix H: 



H -> U- l KU (5.5) 

where A is (block) diagonal and contains the eigenvalues. The Haar measure is then fac- 
torized as follows 



N 

dH = dU J({Xi}) Y[d\i (5.6) 
i=i 

where dU depends only on the eigenvectors and can be integrated out to give a trivial 
constant in front of the integral. J({Xi}) is the Jacobian of the similarity transformation 
O- It is given by j2S| 



MM) ~ II l A * - A / ( 5 - 7 ) 

i<j 

for the Wigner-Dyson ensembles labelled by (3. 

For a gaussian potential V(H) = \H 2 the partition function then takes the form 



JPp({X i })dX 1 ...dX N = Jl[\X i -X j fe-l^^ x Ux 1 ...dX N 



(5. 
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This model is easily solvable. This is the strength of the random matrix description of 
disordered systems. 

The Jacobian in (|5.7|) has the form of a Vandermonde determinant. If we rewrite eq. (15. 7j) 
so that the Jacobian is in the exponent, it gives rise to repulsive eigenvalue correlations in 
the form of a logarithmic pair potential. Such correlations are called geometrical, because 
they arise only from the Jacobian. In the absence of the invariance (|5.4|) . the eigenvalues 
are uncorrelated and follow a Poisson distribution. 



6 The correlators 



The basis for comparing random matrix predictions to experimental or numerical measure- 
ments are the eigenvalue correlation functions. These determine the statistical properties 
of the ensemble. 



The fc-point correlation function pk(Xi, Xk) is defined as 



Pk(Xi, A fe ) = , / J] dXjP({X±, X N }) (6.1) 

1 '' J j=k+l 

Pk(Xi, Xk)dXi...dXk denotes the probability of finding any k eigenvalues in the inter- 
vals dXi around the points Aj (z = l,...,k). The 1-point function is just the density of 
eigenvalues. 

Calculation of the fc-point correlators can be performed exactly by rewriting the Ja- 
cobian as a product of Vandermonde determinants of a set of polynomials orthogonal with 
respect to some measure /(A). This can easily be done just using the properties of deter- 

— /3N A 2 

minants. For instance, for a gaussian matrix model, /(A) = e 57* and the polynomials 
are the Hermite polynomials: 



f +0 ° — /3N A 2 

/ H m (X)H n (X)e dX = h n 5 mn (6.2) 

J — oo 

where h n is a normalization. All the classical orthogonal polynomials appear in this context 
for various choices of the function /(A). 
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The procedure of calculating correlators using orthogonal polynomials was reviewed in 
Mehta's book [2E], arid was a big step forward for random matrix theory at the time. The 
formula for the /c-point correlator can be summarized by the formula 



p /3jfc (Ai, A fe ) = qdet[Qjv/?] (6.3) 

where qdet denotes a quaternion determinant. 4 The matrix elements of Qn/3 are de- 
termined by universal translation invariant kernels j2S] in the large N limit and on the 
unfolded scale. Therefore, eigenvalue correlators are universal and determined only by 
symmetry. To obtain this universal behavior, one has to unfold the spectrum. This means 
rescaling the eigenvalues in such a way that they become dimensionless, by removing the 
dependence of the non-universal density of eigenvalues pi(A). Let's define dimensionless 
variables x« by 



Xi = Xi(Xi) = I pi(X')dX' (6.4) 



The unfolded correlators are obtained by requiring that the differential probabilities should 
be the same in the old and in the new variables: 



Pk(Xi, Xk)dXi...dX k = X k (x!, x k )dx 1 ...dx k (6.5) 
Note that Xi(x) = 1 by construction. 

In practice, however, it is sufficient to do the rescaling in a small region of the spectrum, 
where we are interested in studying the correlations. If this region is centered, say, at the 
origin, we put Zi = Aj/A where A = pi(O) -1 is the average level spacing at the origin 
(note that this quantity is a function of N). This amounts to magnifying the region on the 
scale of the average level spacing, simultaneously with taking the large N limit in which 
the spectrum becomes dense. This so-called microscopic /c-point function (we will call it 
Ps,k{ z ii •••) z k) in accordance with standard usage) is then given by 



ps,k(zi, ~.,Zk) = iv lim A h p k (Azi, Az k ) (6.6) 

N— >oo 



A quaternion determinant of an N x N self-dual quaternion matrix Q is defined as the square root of 
the determinant of the 2N x 2iV matrix obtained by writing each quaternion as a 2 x 2 matrix. 
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Figure 4: Measurements of the spectrum of the Dirac operator in SU(2) lattice gauge theory 
are compared to predictions from random matrix theory. The figure shows the distribution 
of the smallest Dirac eigenvalue P(A m i n ), microscopic spectral density ps{z) (for a 10 4 
lattice) and the microscopic spectral two-point function ps(x,y) (for an 8 4 lattice). The 
histograms represent lattice SU(2) data; the dashed lines are analytical predictions from 
RMT. Reprinted with permission from ref. j22j. Copyright 2005 by the American Physical 
Society. 

where we simultaneously take the limit iV — > oo, keeping the Zi fixed. The new correla- 
tors ps,k(z\, Zk) are independent of the random matrix potential within a large class of 
potentials. The microscopic spectral correlators can be measured in computer simulations 
(see Fig. |3] for an example). 

The cluster functions are relevant for computing spectral observables. A cluster function 
is defined as the connected part of a general /c-point function. (For comparison, for an 
uncorrelated (Poisson) distribution, the connected part vanishes.) For instance, before 
unfolding the two-point function takes the form 

p2,/3(Ai, A 2 ) = pi, / 3(Ai)pi i/3 (A 2 ) - T 2 ,p(\i, A 2 ) (6.7) 

where the 2-level cluster function is given by T 2j/ g(Ai, A 2 ). After unfolding, both the 2- 
point function and the 2-point cluster function depend only on the difference r = x\ — x 2 , 
and we have 

X 2 , p (r) = 1 - Y 2 , p {r) (6.8) 
where X 2> p{r) is the unfolded 2-point function and Y 2j p(r) the unfolded cluster function. 
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7 Spectral observables 



In comparing random matrix theory with experimental and numerical results, we need to 
transform the /c-point functions into statistical spectral observables that can be compared 
to data. 

In this paragraph we closely follow ref. We will consider an energy spectrum that 

is the result of a measurement or of a numerical calculation. The measured values (and 
the random matrix eigenvalues) will be denoted by an ordered sequence {E ly E N }. We 
define a spectral function 



N 

S(E) = J2^E-E n ) (7.1) 

71=1 

and a cumulative spectral function or staircase function 



V (E) = / S(E')dE' = J20(E- E n ) (7.2) 

J n=l 

where 9 is the step function. The cumulative spectral function contains a smooth and a 
fluctuating part: 



V (E) = £(E) + Vfi(E) (7.3) 
where £(-E) is the cumulative mean level density 



£ = = / Pl {E')dE' (7.4) 

J ~ oo 

To obtain the smooth part of the experimentally obtained staircase function, one may fit 
a polynomial to it (see Fig. |SJ). 

The unfolding procedure is identical to the one we performed in RMT in the previous 
paragraph. After unfolding, £ 4 = £(Ei), the staircase function is expressed as 



21 




600 700 800 900 

E [kHz] 




700 702 704 

E [kHz] 



Figure 5: Example of an experimentally obtained staircase function for a spectrum of 1428 
elastomechanical eigenfrequencies of a resonating quartz block. Due to the high number 
of levels, the staircase function appears as a smooth line. The smooth part £(E) is a 
polynomial whose coefficients were found by a fit. The bottom part shows a small section 
of the staircase function. Adapted from Ellegaard et al. [2E] by Guhr, Miiller-Groeling and 
Weidenmiiller j2J. Copyright 2005 by the American Physical Society and with permission 
from Elsevier. 
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(7.5) 



where £ is the dimensionless variable defined in (|7.4|) . Note that the mean level density of 
the unfolded spectrum (i.e. the derivative of the smooth part of the step function) is unity, 
as we have removed the non-universal dependence. 

We will now discuss how a few spectral observables are obtained from the correlation 
functions. 

To study long-range correlations a common statistic is the level number variance £ 2 (L). 
If f] denotes the number of levels in the interval [£, £ + L] in the unfolded spectrum, the 
number variance is defined by 



£ 2 (L) = <rf ) - (r)) 2 = (f ) - L 2 (7.6) 

where the angular brackets (...) denote the average with respect to £. In an interval of 
length L one expects on average L ± a/S 2 (L) levels. 

The number variance is given in terms of correlators by 

H%L) = L — 2 [ (L - r)Y 2tP {r) dr (7.7) 
where Y2${r) is the unfolded 2-level cluster function. 

Another long-range statistic is the spectral rigidity A3. It is defined as the least square 
deviation of the unfolded cumulative spectral function (staircase function) from the best 
fit to a straight line: 

Aa(L) = L-'( (fj(0 - «' + b)f df) (7.8) 

where, like in the number variance, the angular brackets (...) denote the average with 
respect to £. The spectral rigidity can similarly to £ 2 (£) be expressed as an integral 
involving Y^ir). 
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To study fluctuations in the spectrum on a short scale (a few level spacings) we can study 
the nearest neighbor spacing distribution p(s). This is the probability density for two 
neighboring levels £ n and being a distance s apart. The calculation of p(s) is non- 
trivial and involves all correlation functions pk with k > 2. An excellent approximation is 
given by the Wigner surmise: 

pp(s) = a /3 s /3 exp(-& ( gs 2 ) (7.9) 

where ap, bp are /3-dependent constants. Note the level repulsion factor s 13 at small s. See 
Fig. El for an example of measurements of these statistical observables. 

For comparison of measurements to RMT predictions to make sense, we have to make 
the assumption of ergodicity. This means that we assume that the ensemble average in 
the theoretical prediction of RMT is equal to the running average over the sequence of 
measurements on a single sample: 

(f(E)) ens = (f(E)) meas (7.10) 
where f(E) denotes any function of the eigenvalues. 

The observed spectral fluctuations in the systems we have discussed in the introduction 
show, in many instances, an impressive agreement with random matrix theory predictions. 
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Figure 6: Statistical observables in the Dirac spectrum of SU(2) gauge theory. The figures 
show a comparison of lattice data and RMT predictions for the number variance £2, the 
A3 statistic, and the nearest neighbor spacing distribution p(s). The two columns show 
different implementations of fermions on the lattice (left: Wilson fermions, right: Kogut- 
Susskind fermions). GOE, GUE and GSE stand for the Gaussian Orthogonal, Unitary, and 
Symplectic Ensemble, respectively. Reprinted from j2H] with permission from Elsevier. 
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8 Lie groups, algebras, and root lattices 



In this lecture we will first present some preliminary material leading up to the definition of 
symmetric spaces. Assuming that most of the audience is more familiar with this concept 
than the average physicist, we will be as brief as possible. 

As already mentioned, the reason we are interested in symmetric spaces in connection with 
RMT is that random matrix ensembles are identified with symmetric spaces. As we will 
see, symmetric spaces (SS) have well-known properties jUlE] and much of the theory for SS 
can be used in physical problems where RMT is applicable. We will give a few examples 
of such usage in the next lecture. 

We will start by reminding the reader of some basic definitions concerning Lie algebras 
and root spaces. A Lie algebra G is a linear vector space over a field F. Multiplication in 
the Lie algebra is given by the Lie bracket [A, Y]. It has the following properties: 

[1] UX,Ye G, then [A, Y] e G, 

[2] [A, aY + fiZ) = a[X, Y] + {3[X, Z] for a, f3 E F, 

[3] [X,Y]=-[Y,X], 

[4] [A, [Y, Z]] + [Y, [Z, A]] + [Z, [A, Y]] = (the Jacobi identity). 

The algebra G generates a group through the exponential mapping. A general group 
element is 



M = explj2 fX i I 5 t* E. F, A; e G (8.1) 



where ti are parameters (coordinates). We define a mapping adA from the Lie algebra to 
itself by adA : Y — > [A, Y\. The mapping A — > adA is a representation of the Lie algebra 
called the adjoint representation. It is easy to check that it is an automorphism (i.e. that 
it preserves algebraic operations): it follows from the Jacobi identity that [adAj, adA,,] = 
ad[Aj, Xj]. Suppose we choose a basis {A;} for G. Then 



adXiiXj) = [X u Xj] = C%X k (8.2) 

where we sum over k. The C% are real structure constants. The structure constants define 
the matrices M of the adjoint representation through (Mi)jk = C\ k . 
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An ideal I is a subalgebra such that [G,I] C I. A simple Lie algebra has no proper ideal. 
The semisimple algebras are built from the simple ones. In any simple algebra there are 
two kinds of generators. 

(1) There is a maximal abelian subalgebra, called the Cartan subalgebra H = {Hi, H r } 
such that 

[H,H 3 } = (8.3) 

If we represent each element of the Lie algebra by an n x n matrix, then [Hi, H 3 ] = means 
the matrices Hi can all be diagonalized simultaneously. Their eigenvalues /ij are given by 

Hi\n) = tain) (8.4) 

where the eigenvectors are labelled by the weight vectors fi = (pi, /i r ). A positive weight 
is a weight whose first non-zero component is positive. 

(2) There are raising and lowering operators denoted E a such that 

[H h E a ] = a t E a , [E a , E_ a ] = a { H { (8.5) 

Here a is an r-dimensional vector called a root a = (ai, a r ) and r is the rank of the 
algebra. For each root a^, there is another root —a^ and a corresponding eigenoperator 
E- a . The roots form a lattice in the space dual to the Cartan subalgebra. A subset of the 
positive roots span the root lattice. These are the simple roots. Their number is equal to 
r, the rank of the algebra. All the weights of a representation can be obtained by acting 
on the highest weight with lowering operators in all possible ways. 

One can prove the following relation between roots and weights: 

^ = -(p-q) (8-6) 

cr 

where p, q are positive integers such that E a \fi+pa) = 0, E- a \n — qa) = 0, i.e. they define 
the distance of \fi) to the upper and lower end of the ladder. 
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Eq. ()8.6|) implies that the possible angle between two root vectors of a simple Lie algebra 
is limited to multiples of ~ and j. Therefore, there is a finite set of possible root lattices. 
Equation ()8.6j) permits a classification of all complex semisimple algebras. The classical 
Lie algebras SU(n + 1, C), SO(2n + 1, C), Sp(2n, C) and SO(2n, C) correspond to root 
systems A n , B n , C n , and D n , respectively. In addition there are five exceptional algebras, 
but these are not relevant for random matrix theory because they have a finite n. 

The root systems for these four infinite series of classical non-exceptional Lie groups can 
be characterized as follows [HO] (denote the r-dimensional space spanned by the roots by 
V and let {ex, ...e n } be a canonical basis in R n ): 

An-i'. Let V be the hyperplane in R n that passes through the points (1, 0, 0, ...0), (0, 1, 0, 0), 

(0,0, ...,0, 1) (the endpoints of the i — 1, ...,n). Then the root lattice contains the 
vectors {e^ — ej,i ^ j}. 

B n : Let V be R n ; then the roots are {±ej, ie, ± e^, i ^ j}. 
C n : Let V be R n ; then the roots are {±2ei, ±e, ± e,j, i ^ j}. 
D n : Let V be R n ; then the roots are {±ej ± Sj, i ^ j}. 

The root lattice BC n , that we will discuss in conjunction with restricted root systems, is 
the union of B n and C n . It is characterized as follows: 

BC n : Let V be R n ; then the roots are {±e,, ±2e iy ie, ± ej, i ^ j}. 

Roots of length 1, v^2, and 2 are called short, ordinary, and long roots, respectively. Each 
of the complex algebras in general has several real forms associated with it. 5 We will 
define these shortly. Eq. ()8.6|) expresses invariance of the root lattice under reflections in 
the hyperplanes orthogonal to the roots (the Weyl group). If /i is a weight or root, so is 



= /i — a (8.7) 

cr 

The relation (|8.6|) determines the highest weights of all irreducible representations. Setting 
p = 0, choosing a positive integer q, and letting a run through the simple roots, a = a 1 (i = 
1, ...,r), we find the highest weights fi l of all the irreducible representations corresponding 
to the given value of q [SD] • 

5 Also symmetric spaces have real forms, but they will not be discussed here. These are pseudo- 
riemannian symmetric spaces, i.e. they have a non-definite metric. 
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9 Cosets 



In general, a symmetric space can be represented as a coset space of some Lie group G 
with respect to a symmetric subgroup H. The (left) coset space G/H is the set of subsets 
of G of the form gH (g G G): 



G = g H + gi H + ... + g n H (9.1) 

Every element g G G can be written uniquely as g = gfij for some ft £ G and some 
hj G H. The coset can be identified with the set of group operations {g , g n }. The 
coset corresponds to a manifold of dimension dimG — dimif , as we will see in the example 
below. 

Suppose G is represented by matrices acting transitively on a space V (Gv = V for any 
v G V) and Hv = v for some v (then H is called the isotropy subgroup at the point v ). 
Then there is one-to-one correspondence between the elements in V and those in G/H: 
gHv = gv = v. 

Example: The 50(2) subgroup of 50(3) is the isotropy subgroup at the north pole of 
a unit 2-sphere imbedded in 3-dimensional space, since it keeps this point fixed. On the 
other hand, the north pole is mapped onto any point on the surface of the sphere by 
elements of the coset 50(3)/50(2). 

The SU(3) algebra is defined by the commutation relations 



[Li, Lj] — -eij k L k (9.2) 
where 2 e ijk are structure constants. A matrix representation of this algebra is given by 





Li = - 1, L 2 = -\ 00 , L 3 = 77 -1 00 (9.3) 




The subgroup 50(2) is generated by L 3 . This subgroup keeps the north pole fixed: 
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(9.4) 



The remaining group generators define the coset space SO(3)/SO(2). In terms of the real 
coordinates t 1 , t 2 , an element in this coset space takes the form 



M = exp (Elit%) 



WTWf 




;i nv /(ti)2 + (f 2 ) 2 



cos 





2 +(t 2 ) 2 




f(i 2 ) 2 


sin-v/C* 1 ) 


2 +(t 2 ) 2 




f(i 2 ) 2 


Vit 1 ) 2 


+ (t 2 



(9.5) 



The last equation is the equation for a 2-sphere. When the coset space representative M 
acts on the north pole the orbit is exactly the 2-sphere: 




Because of this one-to-one correspondence, the coset space SO(3)/SO(2) can be identified 
with a unit 2-sphere imbedded in 3-dimensional space. 



10 Symmetric spaces 



Suppose G is a compact simple Lie algebra. A linear automorphism o ^ 1 of the Lie 
algebra G onto itself such that a 2 = 1 is called an involutive automorphism or involu- 
tion. This means the eigenvalues of a are ±1, and a splits the algebra G into orthogonal 
eigensubspaces corresponding to these eigenvalues: G = K © P where 
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a(X) =X for X e K, a(X) = -X for X e P 



(10.1) 



K is a subalgebra, but P is not. From eq. (jlO.ip . the following commutation relations hold: 



[K, Kl C K, [K,PlcP, [P,PlcK 



(10.2) 



A subalgebra K satisfying (|l().2jl is called symmetric. If we now multiply the elements in 
P by % (this is called the "Weyl unitary trick"), we construct a new noncompact algebra 
G* = K © iP . This is called a Cartan decomposition of G*, and K is a maximal compact 
subalgebra. The coset spaces G/K and G* / K are symmetric spaces of compact and non- 
compact type, respectively. 

Example: Suppose G = SU(3, C), the group of 3 x 3 unitary complex matrices with unit 
determinant. The algebra of this group consists of eight complex antihermitean traceless 
matrices Xi, i = 1,...,8. Let us take a representation Xi = where Tj denote the 
Gell-Mann matrices known to physicists. 

An involution that splits the SU(3, C) algebra in two subspaces K, P defined as above is 
given by complex conjugation a = K. This involution splits the algebra {Xi, ...X 8 } into 
real and pure imaginary matrices. In the Gell-Mann representation, 



K — {X 2 , X 5 , X 7 } 




{X\,X 3 , X4, Xq, x$} 






2V3 




K is the compact subalgebra SO(3,R) consisting of real, skew-symmetric and traceless 
matrices (this can easily be checked by putting X 2 = L 3 , X 5 = L 2 , X 7 = Li and comparing 
with the 5*0(3) commutation relations [Li, Lj] = ^eijkL k ), and P is the subspace of matrices 
of the form %T , where T is real, symmetric, and traceless. The Cartan subalgebra is given 
by {X 3 ,X 8 }CP. 
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By the Weyl unitary trick we now obtain from G the non-compact algebra G* = K ©z'P, 
where iP is a subspace of real, symmetric, and traceless matrices — T. The entire Lie 
algebra G* consists of 3 x 3 real matrices of zero trace, and generates the linear group of 
transformations SL(3,R). 

The coset space G/K = SU(3,C)/ SO(3, R) is a symmetric space of compact type, and 
the related symmetric space of non-compact type is G*/K = SL(3,R)/SO(3,R). 

Note that the tangent space of G/K (or G*/K) at the origin (identity element) is spanned 
by the subspace P (or iP, respectively) of the algebra. Let's denote by P = e p the 
exponential of any point in the algebra subspace spanned by the set P (such a point is a 
linear combination of the generators in this subspace). When K is a connected subgroup, 
P is isomorphic to G/K. In general P is not a subgroup. However, one can show that if 
p G P, then also p' = kpk" 1 G P. This defines a transitive group action on P. Also, if K 
is compact, every p G P is conjugate with some element in the Cartan subalgebra: 



P 



khk~ 



(10.4) 



This is called spherical decomposition. It defines the angular coordinate k and the spherical 
radial coordinate h of the point p G P. In plain language, every matrix in the coset space 
G/K or G*/K can be diagonalized by a similarity transformation by the subgroup K. 

Example: In the adjoint representation, which has the same dimension as the group, the 
complex symmetric matrices in G/K = SU(3,C)/SO(3) ~ P = e p can be diagonalized 
by the group K = SO (3) to the form 



/ 1 ... 



p = khk~ 



it H 



V 



,it-ct 



-it-a 



Jt-(3 



10.5) 



e-it-7 j 



where we have written the factor of % multiplying the generators in the Cartan subalgebra 
explicitly in the exponent (X 3 = iHi, X 8 = iH 2 ). The vectors ±a, ±/3, ±7 are the 
three pairs of roots of SU(3) (these form a regular hexagon in the plane) and the diagonal 
elements equal to 1 are the exponentials of the zero roots corresponding to the two operators 
in the Cartan subalgebra. 
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11 The metric on a Lie algebra 



A metric tensor can be defined on a Lie algebra. This will be useful for defining the 
curvature of symmetric spaces. Let {X{\ be a basis for the Lie algebra G and let 
denote the structure constants in this basis. The metric tensor on the algebra may be 
defined by the Killing form K(X{, Xj) 

9ij = K(X U X 3 ) = tr(adX l adX J ) = C r is C s jr (11.1) 

The Killing form is symmetric and bilinear. According to Cartan, the Killing form is non- 
degenerate for a semisimple algebra. This means that det^- 7^ 0, so that the inverse of gij, 
denoted by g % \ exists. Since it is also real and symmetric, it can be reduced to canonical 
diagonal form g^ = diag(— 1, —1, 1, 1). 

According to a theorem by Weyl, a simple Lie group G is compact, if and only if the Killing 
form on G is negative definite. Otherwise it is non-compact. 

Example: We have already written down the commutation relations of the compact group 
SO(3,R) in the form [Li,Lj] = CfjLk = \eijkLk- We can renormalize the generators so 
that the entries of the metric are unity. The commutation relations then take the form 

[Li, L 2 ] = -±L 3 , [L 2 , L 3 ] =-^L u [L 3 , L x } = -^L 2 (11.2) 

We can read off the structure constants and then, using eq. (|ll.lj) . compute the components 
of the Killing form. In this normalization it is 




(11.3) 



(In the unrenormalized form it is g^ = —§<%.) This is the metric in the 5*0(3) algebra. It 
is negative, because the group is compact. 

The generators of the non-compact group SO (2, 1; R) obey the commutation relations 

[E l5 S 2 ] = - 1 E 3 , [S 2 , S 3 ] = --tjSi, [S 3 , Si] = -^S 2 (11.4) 
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In the same way, using again eq. (Jll.lj) we compute the matrix elements gij. The result is 




(11.5) 



we have labelled the rows and columns in the order 3,1,2. The generator S 2 makes up the 
compact subalgebra K. 



12 The metric on a symmetric space 

The definition of the metric can be extended to an arbitrary point of a symmetric space. 
At the origin, the metric is defined by restricting the metric in the algebra to the tangent 
space (remember that the latter is spanned by the generators in the subspace P or iP of 
the whole algebra G or G*). 

Example: The metric in the subspace P of SO(3,R) is obtained by excluding the row 
and column corresponding to the generator in K, keeping the ones in P: 




(12.1) 



Similarly, the metric in the subspace iP of SO(2, 1; R) is obtained by excluding the row 
and column corresponding to the compact generator in K and keeping the ones in iP: 




(12.2) 



Since the group acts transitively on the symmetric space, we can then use a group trans- 
formation to map the metric to an arbitrary point of the SS, using the invariance of the 
line element in local coordinates given by ds 2 = gijdx % dxK 

Note that if a positive metric is required on a compact symmetric space, we can use minus 
the Killing form, which sometimes is more natural. 

Example: The line element ds 2 on the radius-1 2-sphere isomorphic to the symmetric 
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space SO(3, R)/S0(2) in polar coordinates is ds 2 = d9 2 + sin 2 ^^ 2 . The metric at the 
point (9, 0) is 

= ( J si „\ ) • <r = ( I J*g ) < 12 - 3 > 

where the rows and columns are labelled in the order 9, <ft. 

The line element ds 2 on the hyperboloid SO(2, 1; R)/SO(2) in polar coordinates is ds 2 = 
d9 2 + sinh 2 # dcj) 2 . The metric at the point (6, (ft) is 

= ( I sil 2 9 ) ' = ( J sinh% ) ^ 12 - 4 ) 

13 Real forms and the metric 

The form of the metric depends on the basis of the algebra. A complex Lie algebra G c is 
given by 

G c = ciR i + c ° Ea ( c *> ° a com P lex ) t 13 - 1 ) 

i a 

where Ho = {Hi} is the Cartan subalgebra and {E± a } are the pairs of raising and lowering 
operators. A real form of the same algebra is obtained by taking the coordinates c 1 , c a to 
be real numbers, i.e. 

G = ^2dHi + ^2 c a E a (c\ c a real) (13.2) 

i a 

We can choose different basis vectors in this real algebra. The form of the metric will 
change accordingly. 

The metric corresponding to the basis {Hi, ±E a } is not diagonal. It has the form 
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/ 1 



9ij 



1 

1 



\ 



1 

1 0/ 



;i3.3) 



By recombining the basis vectors into 



K 



(E a — E_a) 

7i 



(Eg + E_ a ) 

V2 



13.4) 



the metric takes the diagonal form 



/ 1 



9ij 



\ 



1 
-1 



1 





\ 



■1 / 



;i3.5) 



This is called the normal real form with a diagonal metric. Here the labelling of the rows 
and columns is such that entries with a minus sign correspond to the generators of the 
compact subalgebra K, the first r entries equal to +1 correspond to the Cartan subalgebra, 
and the remaining ones to the operators in iP not in the Cartan subalgebra. 

The compact real form is obtained from the normal real form by the Weyl unitary trick in 
reverse. That is, we choose the basis to be 



K 



(Eg — E-g) 

7! 



iHi 



i{Eg + E_ 

V2 



;i3.6) 
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The metric tensor is then 



/ -1 



9ij 



\ 



;i3.7) 



Example: We will use as an example the well-known SU(2, C) algebra with Cartan 
subalgebra Hq — {J3} and raising and lowering operators {J±}- 



2V2 

where 



^3 = ^3, J ± = i(n±zr 2 ) (13* 



We have chosen the normalization such that the non-zero entries of are all equal to 1: 



[Js, J±] = ±^ J±, [J + ,J-] = j= 2 J 3 (13.10) 

From the commutation relations we can read off the structure constants and determine the 
corresponding metric. It is non-diagonal: 



/ 1 
^•=001 
\ 1 

where the rows and columns are labelled by 3, +, — respectively. To pass now to a diagonal 
metric, we use the recipe in eq. (I13.4|) 



13.11) 
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£i = J -^t = ( 13 - 12 ) 

V — J+-J- _ i _ 



The commutation relations then become 



[E 1; E 2 ] = -^E 3 , [E 2 , E 3 ] = -^Ei, [E 3 , Ex] = ^E 2 (13.13) 

These commutation relations characterize the algebra SO(2, 1; R). From here we find the 
non-zero structure constants Cf 2 = — Cfi = C 23 = — C\ 2 = — Cfi = C*i 3 = — ^75 and the 
diagonal metric of the normal real form with rows and columns labelled 3, 1, 2 (in order to 
comply with the notation in eq. (|13.5[) ) is 



/ 1 
9H = 1 
\ -1 

which is to be compared with eq. ()13.5j) . According to eq. (|13.4j) . the Cartan decomposition 
of G* is G* = K©zP where K = {S2} and iP = {E3, Ei}. The Cartan subalgebra consists 
of E 3 . 

Finally, we arrive at the compact real form by multiplying E 3 and Ei with i. Setting 
fEi = Ei, E 2 = S 2 , «E 3 = E 3 the commutation relations become those of the special 
orthogonal group: 



(13.14) 



[E 1 ,E 2 ] = -^E 3 , [E 2 ,E 3 ]=-iE!, [E 3 ,E 1 ] = -^E 2 (13.15) 

The last commutation relation in eq. (J11.4j) has changed sign whereas the others are un- 
changed. C| x , Cf 3 , and consequently g 33 and g±\ change sign and we get the metric for 
SO(3,R): 



9ij 
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(13.16) 



This is a compact metric and SO(3,R) is the compact real form of the complex algebra. 
The subspaces of the compact algebra G = K © P are K = {E 2 } and P = {E 3 , Ei}. 

To summarize, real forms are obtained by using different combinations of basis vectors 
(generators) and by using the Weyl unitary trick. This changes the commutation relations 
and the form of the metric. In general, a semisimple complex algebra has several distinct 
real forms: one compact and several non-compact ones distinguished by their character 
(trace of the metric) . 



14 Obtaining all the real forms of a complex algebra 



In the previous section we saw how to construct the compact real form of an algebra. 
To classify all the real forms of any complex Lie algebra, it suffices to enumerate all the 
involutive automorphisms of its compact real form. If G is the compact real form of 
a complex semisimple Lie algebra G c , G* runs through all its associated non-compact 
real forms G*, G'*, ... with corresponding maximal compact subgroups K, K', ... and 
complementary subspaces iP, iP', ... as a runs through all the involutive automorphisms 
of G (W denotes the Weyl trick) 



ai+w G* = K © iP 
/ 

G c -> G n ±F G'* = K' © iP' (14.1) 

CT3+W 

\ Q//* = K // ip // 

One can show [31] (Ch. VII), that it suffices to consider the following three possibilities (or 
combinations thereof) for a: a\ — K (complex conjugation), cr 2 = I P ,q and 03 = J PiP where 




(14.2) 



and I p denotes the p x p unit matrix. 

An arbitrary involution crj mixes the subspaces K and P of the compact real form and 
splits the algebra in a different way into K' © P'. The non-compact real forms are then 
obtained through the Weyl unitary trick. 
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Example: The algebra SO(3,R), belonging to the root lattice B\ is spanned by the 
generators L±, L 2 , L 3 given in sectional A general element of the algebra is 




(14.3) 



This splitting of the algebra is caused by the involution 7 2 ,i acting on the representation: 



1 



t 3 t 2 



-t 3 t 1 



1 

2 



t 3 -t 2 
-t 3 -t 1 
t 2 t 1 



;i4.4) 



and it splits it into SO (3) = K © P = SO (2) © SO (3) /SO (2). Exponentiating, the coset 
representative is a point on the 2-sphere 



M 



x 

y I ; x 2 + y 2 + z 2 = l 



(14.5) 



By the Weyl unitary trick we now get the non-compact real form G* = K©zP: SO(2, 1) 
SO(2) © SO(2, l)/SO(2). This algebra is represented by 



-t a 



t 3 It 2 \ 




1 e * ^ 




It 1 


\-\ 




|e( 


-It 1 ) 









it 2 
it 1 



-it 2 -it 1 



(14.6) 



and after exponentiation of the coset generators 



M 



IX 

iy 

z 



(ix) 2 + (iy) 2 + z 2 = 1 



(14.7) 



The surface in R 3 consisting of points (x, y, z) satisfying this equation is the hyperboloid 
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H 2 . Similarly, we get the isomorphic space SO(l, 2)/SO(2) by applying Ji 2 : SO(l,2) = 
K © iP = SO(2) © SO(l, 2)/SO(2) and in terms of the algebra 




15 The classification of symmetric spaces 

The reason we have discussed how to obtain all the real forms of a complex algebra is 
that we want to understand how to obtain all the riemannian symmetric spaces associated 
with the simple Lie groups. These are namely exactly the symmetric spaces appearing 
as integration manifolds in hermitean random matrix theory. (There are also pseudo- 
riemannian symmetric spaces, but these are not too interesting here.) Once we have the 
real forms, the symmetric spaces are defined as the spaces corresponding to the exponential 
mapping of the subspaces P and iP of the various real forms of the complex algebra. In 
addition, if G is a compact group, G and G c /G are also such riemannian symmetric spaces. 

To summarize, the interesting manifolds are obtained by 



• obtaining the compact real form of a simple Lie algebra by combining ladder operators 
appropriately; 

• operating with all possible involutions on the resulting algebra, thereby obtaining all 
the possible symmetric subalgebras; 

• forming the corresponding non-compact algebras by the Weyl unitary trick; 

• then forming pairs of symmetric coset spaces G/K, G*/K, G'/K', G'*/K',... corre- 
sponding to the various compact symmetric subgroups and in addition, including the 
symmetric spaces G and G c /G. 



16 Curvature 



A curvature tensor with components K t j kl can be defined on the manifold G/K or G*/K. 
It is given by 
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R-ijkXn — [Xi, [Xj,Xk]] — c™ m cyix n (16.1) 

where {Xi} is a basis for the Lie algebra. If {X, Y} is an orthonormal basis for a two- 
dimensional subspace S of the tangent space at a point p (assumed to have dimension > 2), 
the sectional curvature at the point p is 

K = g{[[X,Y],X] % Y) (16.2) 

For a two-dimensional manifold, this is just the Gaussian curvature. If the manifold has 
dimension > 2, (J16.2J1 gives the sectional curvature along the section S. These equations, 
together with the commutation relations for K and P, show that the curvature of the 
spaces G/K and G* / K has a definite and opposite sign. To the same subgroup K there 
corresponds a positive curvature space P ~ G/K and a dual negative curvature space 
P* ~ G*/K. 

Example: We can use the example of SU(2) to see that the sectional curvature is the 
opposite for the two spaces G/K and G*/K. If we take {X, Y} = as the basis 

in the space iP and {S 3 , Ei} (Ej = zEj) as the basis in the space P, we see by comparing 
the signs of the entries of the metrics we computed in eqs. (|13.14|) and (|13.16j) that the 
sectional curvature K at the origin has the opposite sign for the two spaces 5*0(2, 1)/ S0{2) 
and SO(3)/SO(2). 

In addition, there is also a zero-curvature symmetric space X° = G°/K related to X + = 
G/K and X~ = G* / K. The space X° can be identified with the subspace P of the algebra. 
The group G° is a semidirect product of the subgroup K and the invariant subspace P of 
the algebra, and its elements g = (k, a) act on the elements of X° in the following way: 

g(x) = kx + a, k G K, x, a G X° (16.3) 

if the x's are vectors, and 

g(x) = kxk~ l + a, k e K, x,aeX° (16.4) 

if they are matrices. The elements of the algebra P now define an abelian additive group, 
and X° is a vector space with euclidean geometry. In the above scenario, the subspace P 
contains only the operators of the Cartan subalgebra and no others: P = H , so that both 
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K and P in this case are subalgebras of G°. The algebra G° = K © P is non-semisimple 
because the subgroup P is an abelian ideal ([P, P] = and [K, P] C P). Note that K and 
P still satisfy the commutation relations in eq. ()10.2|) . By this equation, R^ k = for all 
the elements X G P. 

Even though the Killing form on non-semisimple algebras is degenerate, it is trivial to find 
a non-degenerate metric on the symmetric space X°. 

Example: An example of a flat symmetric space is E 2 j SO(2), where G° = E 2 is the 
euclidean group of motions of the plane R 2 : g{x) = kx + a, g = (k,a) G G° where 
k G K = SO (2) and a G R 2 . The generators of this group are translations Pi, P 2 G H = P 
and a rotation J G K satisfying 

[P u P 2 ] = 0, [J, Pi] = -e^Pj, [J, J] = (16.5) 

in agreement with eq. ()10.2|) defining a symmetric subgroup. The abelian algebra of trans- 
lations ^21=1 t % Pii t l ^ R-j is isomorphic to the plane R 2 , and can be identified with it. 

The Killing form for E 2 is degenerate: 




(16.6) 



Therefore we don't take the Killing form as the metric on E 2 j SO(2). Instead, we can use 
the Euclidean metric 

J) (16.7) 

on the entire symmetric space. 

We remark that the zero-curvature symmetric spaces correspond to the integration mani- 
folds of many known matrix models with physical applications. 
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17 Restricted root systems 



The restricted root systems play an important role in connection with matrix models 
and integrable Calogero-Sutherland models. Here we will only describe very briefly how 
restricted root systems are obtained and how they are related to a given symmetric space. 
Due to lack of space, we will not give any example of the construction of such a root space. 
A concrete example was worked out in subsection 5.2 of reference 

Real forms of a complex algebra share the same root system with the latter. This is because 
they correspond to the same set of raising and lowering operators and the same Cartan 
subalgebra. One can also associate a root system to a symmetric space G/K. If the root 
system of the group G has rank r, the rank of this restricted root system may be different, 
say r' . 

Example: The algebra SU(p, q; C) (p + q = n) is a non-compact real form of SU(n, C). 
They share the same rank-(n — 1) root system A n _i. The restricted root system of the 
symmetric space SU(p,q;C)/(SU(p) Cg> SU(q) ® U(l)) is BC r /, where r' = min(p,q). 

In general the restricted root system will be different from the original, inherited root 
system if the Cartan subalgebra is a subset of K. The procedure to find the restricted root 
system is then to redefine the Cartan subalgebra so that it lies partly or entirely in P (or 
in iP, if appropriate). This is possible if we 

• find a new representation of the original Cartan subalgebra Ho corresponding to the 
original root lattice. This corresponds to a Weyl reflection of the root lattice and 
can be achieved by a permutation of the root vectors. In practice it amounts to 
permuting the diagonal elements of the original H^s. 

• do this in such a way that a maximal number r' of commuting generators are in the 
subspace P. The new Cartan subalgebra A has the same number of generators as 
Ho (this number equals the rank of the algebra), but r' of its elements lie in the 
subspace P. r' is called the rank of the symmetric space G/K. 

The new root system is defined with respect to the part of the maximal abelian subalgebra 
that lies in P. Therefore its rank can be smaller than the rank of the root system inherited 
from the complex extension algebra. We can define raising and lowering operators E' a in 
the whole algebra G that satisfy 

[X' t , E' a ] = a[E' a (i;eA nP) (17.1) 
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The roots a[ define the restricted root system. In addition to the sign of the curvature 
of the symmetric space, it is the restricted roots (and their multiplicities) that define the 
Jacobian in the transformation to radial coordinates on the symmetric space. This Jacobian 
is exactly the one we encounter in random matrix theory too, when we diagonalize the 
ensemble of random matrices to obtain a partition function (and correlators) expressed only 
as a function of random matrix eigenvalues. The latter are exactly the radial coordinates. 
In RMT we integrate out the degrees of freedom corresponding to the symmetric subgroup. 

We have discussed the procedure for obtaining the irreducible riemannian symmetric spaces 
originating in simple Lie groups. They are listed in Table El The classification due to E. 
Cartan is based on the root systems inherited from the complex extension algebra. We 
have also explained how the restricted root lattice is defined for each symmetric space. For 
each compact symmetric subgroup K there is a triplet of symmetric spaces corresponding 
to positive, zero, and negative curvature. The zero curvature spaces are isomorphic to 
algebra subspaces P (which are the tangent spaces of the symmetric spaces of positive 
curvature) and are not listed. The root multiplicities pertain to the restricted root systems 
of the pairs of dual symmetric spaces with positive and negative curvature. 

As explained, the real forms of the simple Lie groups do not include all the possible 
riemannian symmetric spaces. The compact Lie group G is itself such a space, and so is 
its dual G c /G (here the algebra G c = G* © iG* is the complex extension of all the real 
forms G*). These are also listed in Table El 



18 Invariant operators on symmetric spaces 

Let G be a (semi)simple rank-r Lie algebra. A Casimir operator Ck (k = 1, r) associated 
with the algebra G is a homogeneous polynomial operator that satisfies 



for all Xi e G. The simplest (quadratic) Casimir operator associated to the adjoint 
representation of the algebra G is given by 



[C k ,Xi] = 




C = gVXiXj 
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Table 2: The classification of irreducible symmetric spaces of positive and negative curva- 
ture originating in simple Lie groups. 



Inherited 
root space 


Restricted 
root space 


Car tan 
class 


G/K (G) 


G*/K (G c /G) 


m Q 


mi 


m s 


A 

Azv-i 


A 

Azv-i 
A N -i 
A N -! 

D^Jq (P>q) 
C q (p=q) 


A 

AI 
All 

AIII 


QTT 1 AT\ 
oU {1\ ) 

SU(N) 

SO(N) 

SU(2N) 

USp(2N) 

SU{p+q) 


SL{N,C) 

SU(N) 
SL(N,R) 

SO(N) 
SU*(2N) 
USp(2N) 

SU(p,q) 


o 
Z 

1 

4 
2 


U 



1 


U 



2(p-q) 


SU{p)xSU(q)xU(l) 


SU(p)xSU{q)xU(l) 




dn 


~D 
D 


oUyZlv + i) 


SO(2N+l,C) 
SO(2N+l) 




n 
U 


o 
Z 


Cn 


Cn 
Cn 

BCq (p>q) 
C q (p=q) 


c 

CI 
CII 


USp(2N) 

USp(2N) 


Sp(2N,C) 
USp(2N) 
Sp(2N,R) 


2 
1 

4 


2 
1 

3 





4(p - q) 


SU(N)xU(l) 
USp(2p+2q) 


SU(N)xU(l) 
USp{2p,2q) 


USp{2p)xUSp(2q) 


U Sp(2p)xU Sp(2q) 


D n 


D N 
Cn 
BCn 


D 

Dill-even 
Dill-odd 


S0(2N) 

SO(W) 


SO(2N,C) 
SO(2N) 
SO*(4N) 


2 
4 
4 



1 
1 




4 


SU(2N)xU(l) 
SO(4N+2) 


SU(2N)xU(l) 
SO* {AN +2) 


SU(2N+l)xU(l) 


SU(2N+l)xU(l) 


Bn {p+q=2N+l) 
D N (p+q=2N) 


B q (p>q) 
D q (p=q) 


BDI 


SO(p+q) 


SO(p,q) 


1 





p-q 


SO(p)xSO(q) 


SO(p)xSO(q) 



where g lJ is the inverse of the metric tensor 6 defined in (111.1)1 and the generators Xi are in 
the adjoint representation (it can be defined in a similar way for any other representation 
of G). 

In general, the product XY makes no sense in the algebra G. The Casimir operators lie 
in the enveloping algebra obtained by embedding G in the associative algebra defined by 
the relations 

X(YZ) = (XY)Z [X,Y]=XY-YX (18.3) 

The number of functionally independent Casimir operators is equal to the rank r of the 
group. 

6 Note that Casimir operators are defined for semisimple algebras, where the Killing form is non- 
degenerate. This does not prevent one from finding operators that commute with all the generators of 
non-semisimplc algebras. For example, for the euclidean group £3 of rotations { Ji, J2, J3} and translations 
{Pi, P2, P3}, P 2 = PiPi an d P J = PiJi commute with all the generators. Also the operators that 
commute with all the generators of a non-semisimple algebra are often referred to as Casimir operators. 
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All the independent Casimir operators of the algebra G can be obtained by making the 
substitution t l —* X,i in the functionally independent coefficients (fik{t l ) of the secular 
equation: 



det (£S G iVM - xin) = ELo(-A) n -V,(i 



Ci(Xi 



18.4) 



where p(Xi) is some representation of the algebra and (fk{t l ) are functions of the real 
coordinates t l . In general there will be r functionally independent coefficients, and r 
functionally independent Casimir operators. 

Example: The secular equation for the SO (3) rank-1 algebra is 



det ft ■ L - \h 



-A t 3 /2 t 2 /2 
-t 3 /2 -A t x /2 
-t 2 /2 -t l /2 -A 



-A) 3 + (-A)V = 



;i8.5) 



As expected, this equation has one functionally independent coefficient, <^i(t) 
only Casimir operator is the square of the angular momentum operator: 



it 2 . The 



C\ ~ L 2 = L\ + h\ + L 2 



18.6) 



obtained by the substitution t l — > Li in (pi(t). The Casimir operator can also be obtained 
from eq. ()18.2|) by using the metric gij = —\5ij for SO (3) given in a previous example. We 
already know from quantum mechanics that 



[L 2 ,L!] = [L 2 ,L 2 ] = [L 2 ,L 3 ] =0 



;is.7) 



The Casimir operators can be expressed as differential operators in the local coordinates 
on the symmetric space: 



_d_ 

dx c 
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where x a are local coordinates [SJ|32] (for example, L x = (r x p) x = —i(yd z — zd y )). 

Expressed in local coordinates as differential operators, the Casimirs are called Laplace 
operators. In analogy with the Laplacian in R n , 



(18 ' 9) 

i=i 

which is is invariant under the group E n of rigid motions (isometries) of R n , the Laplace 
operators on (pseudo-)riemannian manifolds are invariant under the group of isometries 
of the manifold. The isometry group of the symmetric space P ~ G/K is G, since G acts 
transitively on this space and preserves the metric. The number of independent Laplace 
operators on a riemannian symmetric coset space is equal to the rank of the space. 

The Laplace-Beltrami operator on a symmetric space is the special second order Laplace 
operator. It can be expressed as 



^=-^iz^w l > 9 = det9 « (1810) 

Example: Let's calculate the Laplace-Beltrami operator on the symmetric space 50(3) / 50(2) 
in polar coordinates using (|18.1U|) and the metric at the point (6, 4>) 



* = I J 8in»» ) ■ o" = ( I J-h I 



Substituting in the formula and computing derivatives we obtain the Laplace-Beltrami 
operator on the sphere of radius 1: 



A B = d] + co^ d e + sin" 2 ^ d\ (18.12) 

Of course this operator is proportional to L 2 . We can check this by computing L x = 
—i(yd z — zd y ), L y = —i(zd x — xd z ), and L z = —i(xd y — yd x ) in spherical coordinates (setting 
r = 1) and then forming the operator L 2 X + L 2 + L 2 Z , remembering that all the operators 
have to act also on anything coming after the expression for each L\. We find that L 2 in 
spherical coordinates, expressed as a differential operator, is exactly the Laplace-Beltrami 
operator. 
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As we have seen in eq. fjlU.4|) . radial coordinates can be defined on the SS. The adjoint 
representation of a general element H in the maximal abelian subalgebra H' C P follows 
from a form similar to eq. ()10.5)) (with or without a factor of i depending on whether we 
have a compact or non-compact space), but now the roots are in the restricted root lattice. 
For a non-compact space of type P* 



logh = H = q • H 



q • a 



\ 



18.13) 



-q -7} J 



We define q a = q • a. These are the radial coordinates on the symmetric space. 

Example: The rank of the symmetric space SU(2,C)/SO(2) is 1 and the restricted root 
lattice is A\. Absorbing a factor of \[2 (the length of the ordinary roots) into the coordinate, 
the above equation takes the form 



H = 6H 1 = 6 



h = e WH ' 
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18.14) 



The radial coordinate is q = (q 1 ) = 9. 

Example: From Tabled on the symmetric negative curvature space SO(2N, C)/SO(2N), 
the restricted root lattice is of type and the roots are {±ej ± ej,i ^ j}. The q a are 
±qi ± qj- 

In general, a Laplace-Beltrami operator can be split into a radial part A' B and a transversal 
part. The radial part acts on geodesies orthogonal to some submanifold S, typically a 
sphere centered at the origin jjj. 

Example: The radial part of the Laplace-Beltrami operator for the coset space SO(3) / SO (2) 
given in ([18.12)1 is 



A' R = dj + cot6d B 



18.15) 
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The radial part of the Laplace-Beltrami operator on a symmetric space has the general 
form 



£ = J-X^-SLjQ)-SL (j = 0,-+) (18.16) 
JO) Qq a dq a ^ ' V ' 

a=l 

where r' is the rank of the symmetric space, is the Jacobian of the transformation to 
radial coordinates on the SS (to be given below) and m a is the multiplicity of the restricted 
root a. (The multiplicities m a were listed in Tabled) The sum in (J18.16)) goes over the 
labels of the independent radial coordinates q = \ogh(x) = (q 1 , .■■,q r ) where h(x) is the 
exponential map of an element in the Cartan subalgebra. 

The Jacobian in ()18.16|) is given by 



J(°)(g) 


= IW(? a ) ma 


J(~\q) 


= n ae K+(sinh(<f)r 


jC+)(g) 





18.17) 



for the various types of symmetric spaces with zero, negative and positive curvature, re- 
spectively (see [7], Ch. I, par. 5). = y/\g\ where g is the metric tensor at an arbitrary 
point of the symmetric space. In these equations the products denoted IlaeiJ+ are over a ^ 
the positive roots of the restricted root lattice. 7 . 

Example: If the restricted root lattice is of type An with only ordinary roots ej — ej 
(i ^ j), the Jacobian of the zero curvature space is 



^ 0) ({ft}) = nifc-<&r° ( i8 - 18 ) 

i<j 

The absolute value corresponds to a certain choice of Weyl chamber (ordering of the qi) 

m- 

7 Strictly speaking, in the cuclidean case we have not defined any restricted root lattice. The formula for 
the Jacobian (q) for the zero-curvature space is understood as the infinitesimal version of the formula 
pertaining to the negative-curvature space. 
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Example: If the restricted root lattice is of type Cn with long and ordinary roots, the 
positive roots are {e^ ± ej, 2e;}. The Jacobian of the negative curvature space is then 

= Y[ I sinh 2 ^ - smh 2 qj \ m ° Y[smh mi {2q k ) (18.19) 

i<j k 

where we have used sinh(gj — gj)sinh(gj + qj) = sinh 2 gj — sinh 2 g,,-. (The root multiplicities 
were listed in Tabled) 

Since the Laplace operators form a commutative algebra, they have common eigenfunc- 
tions. The eigenfunctions of the radial part of the Laplace-Beltrami operator on a sym- 
metric space are called zonal spherical functions. They play an important role in math- 
ematics as bases for square-integrable functions, not to mention their role as irreducible 
representation functions in quantum mechanics. In the present context we will see that 
they determine the solution of some physical problems where the relevant operator can be 
mapped onto the Laplace-Beltrami operator. 

Suppose the smooth complex-valued function <fi\(x) is an eigenfunction of some invariant 
differential operator on the symmetric space G/K: 

A fe A (x) = 7Afe (A)0 A (x) (18.20) 

The function 4>\(x) is called spherical if it satisfies cf)\{kxk') = 4>\{x) (x G G/K, k G K) and 
if (f)\{e) = 1 (e =identity element). Because of the bi-invariance under K, these functions 
depend only on the radial coordinates h: 

4>\(x) = Mh) (18.21) 

Example: We know from quantum mechanics that the eigenfunctions of the Laplace 
operator L 2 on G/K = SO(3)/SO(2) are the associated Legendre polynomials P;(cos#). 
Setting L x = —i(yd z — zd y ) etc., 

L 2 P,(cos#) = 1(1 + l)P,(cos0) (18.22) 

where cos# is the ^-coordinate of the point P = (x, y, z) on the sphere of radius 1 (in 
spherical coordinates, P = (sin^ cos<f>, sin6 sin0, cos6*)). As we can see, the eigenfunctions 
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are functions of the radial coordinate 9 only. The subgroup that keeps the north pole fixed 
is K = SO{2) and its algebra contains the operator L z = 8$. Indeed, P;(cos#) is unchanged 
if the point P is rotated around the z-axis. 

Following reference [Oj, we introduce a parameter a into the the Jacobians (jl8.17j) for the 
symmetric spaces, 



jW(q) 


= riaeR+' 






J {+) (q) 





18.23) 



The parameter a corresponds to a radius. For example, for the sphere SO(3)/SO(2) it is 
the radius of the 2-sphere. 

The various spherical functions corresponding to the spaces of positive, negative and zero 
curvature are then related to each other by the simple transformations [H| 



^ o) (g)=lim a ^o0i _) (g) 
4 +) (?) = 0i" ) (?)l™ 



18.24) 



and their eigenvalues are given by 



Arf } = -A 2 0i O) 
(-) _ _ „ 



= - p 2 w ( i8 - 25 ) 

AU +) = (-§ + P 2 )0i +) 



where p is the function defined by 



2 
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There is an extensive theory relating to such eigenfunctions [7] , but we will not be able to 
discuss it here. Some important results and references were listed in jj. 



19 A new classification of RMT 

In this third lecture we will discuss three applications of the theory of symmetric spaces in 
random matrix theory. We have not been able to dwell upon the details of all the various 
random matrix ensembles. However, they all have some common features that we have 
already discussed in the introductory part of these lectures: 

• the ensembles are determined by physical symmetries and labelled by a Dyson index 
(3 counting the degrees of freedom of the matrix elements (there is also a boundary 
index a characterising the ensembles that do not have translationally symmetric 
eigenvalue distributions); 

• the probability P(M)dM is invariant under some similarity transformation of the 
matrix M; 

• the random matrix integral can be expressed as a function of random matrix eigen- 
values by diagonalizing the ensemble; 

• the Jacobian due to this diagonalization determines the geometric eigenvalue repul- 
sion characteristic of RMT; 

• one can identify ensembles of random matrices with symmetric spaces. 

We will now discuss the correspondence between matrix ensembles and symmetric spaces 
in more detail through a few examples. 

Example: The circular random matrix ensembles are ensembles of random scattering 
matrices S relating the incoming and the outgoing wave amplitudes in a scattering problem. 
Scattering processes are important both in mesoscopic physics and in many-body problems 
in nuclear and atomic physics. The scattering system is idealized as incoming and outgoing 
scattering channels in which propagation is free, and a compact interaction region where the 
scattering takes place. Let us for simplicity consider scattering in a mesoscopic disordered 
system connected to electron reservoirs through leads. In Fig. a schematic wire of length 
L and width W is shown (the region II is disordered). The wave functions of incoming and 
outgoing electrons in the left and right leads are denoted ^% L (r)- A wave function \E f± (f^ 
can be decomposed as 



53 



X 



I II III 




Figure 7: A schematic picture of a quantum wire consisting of a disordered region (II) of 
length L and width W connected to electron reservoirs through leads (I, III). In each lead 
there are N incoming and N outgoing scattering channels. 



*±(f) = <S> n (r t )e ±ik " x (19.1) 

where <& n (r t ) is the transverse wave function and the integer n = 1,...,N labels the N 
propagating modes or scattering channels. The coordinate x is along the wire. 

The scattering matrix S relates the incoming and the outgoing wave amplitudes of the 
electrons. If we have N propagating modes at the Fermi level, we can describe them by a 
vector of length 2N of incident modes I, V and a similar vector of outgoing modes O, O' 
in each lead, where unprimed letters denote the modes in the left lead and primed letters 
the modes in the right lead. Then the scattering matrix is defined by 




(19.2) 



Flux conservation (|/| 2 + \I'\ 2 = \0\ 2 + \0'\ 2 ), implies that S is unitary. The symmetry 
classes we discussed in connection with Wigner's hamiltonian ensembles are reflected here 
in the circular orthogonal ensemble {(3 = 1, S is unitary and symmetric), the cirular 
unitary ensemble {(3 — 2, S is unitary) and the circular symplectic ensemble ((3 = 4, S 
is unitary and self-dual). (Note that the hermitean Hamiltonians can be related to the 
unitary scattering matrix by S = e lH .) 

Let's look more closely at the circular orthogonal ensemble. Every symmetric unitary 
matrix S can be written as 
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S = U T U 



(19.3) 



where U is a generic unitary matrix. However, this mapping is not one-to-one. If we 
assume that S = U T U = V T V, then it is easy to see that the matrix relating the two 
expressions, R = VU~ l , is unitary and satisfies R T R =10- Hence R must be real and 
orthogonal. Thus we see that the manifold of the unitary symmetric matrices is actually 
the coset U(N)/0(N), due to the above mentioned degeneracy. From the point of view 
of the physical properties of the ensemble nothing changes if we perform the restriction 
to an irreducible symmetric space. 8 Then the manifold becomes SU(N)/ SO(N). The 
integration manifold in the unitary ensemble is simply the group SU(N) without further 
constraint and in the symplectic ensemble, in an analogous manner to the orthogonal 
case, we realize that the manifold coincides with SU(2N)/ Sp(2N). Comparing now with 
Table we see that the integration manifolds of the three circular ensembles are exactly 
the first three coset spaces (described in the Cartan notation as A, AI and All) of positive 
curvature in the list of possible irreducible symmetric spaces. 

Example: As implied in the discussion of hamiltonian ensembles in section Pp(H) and 
the integration measure dH are separately invariant under the transformation 

H -> UHU~\ (19.5) 

where U is an orthogonal, unitary or symplectic N x N matrix, depending on the value of 
(3. It can be shown J2E] that the form of Pp(H) is automatically restricted to the form 

Pp{H) = exp(-atrF 2 + btiH + c) (19.6) 

(a > 0) if one postulates statistical independence of the matrix elements Hij. Note that 
Pp(H) can be cast in the form 

Pp{H) ~ e~ atrH2 (19.7) 

8 In the partition function 

Z~ / dSPp{S) (19.4) 
Jg/k 

extracting such a U(l) factor from the integration manifold just amounts to redefining Z by a constant. 
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by simply completing the square in the exponent. As we know, this is a good choice for 
the RMT potential, because it makes the theory easy to solve. 

However, it is important to notice that the symmetry group of Pp(H) dH is larger and 
consists of rotations by the matrix U like in eq. (j!9.5jl . and addition by square hermitean 
matrices: 



This latter equation tells us that the ensemble is translation invariant. Limiting the dis- 
cussion again to irreducible symmetric spaces, for /3 = 1 we are then dealing with the 
set of real, symmetric and traceless matrices. But this is exactly the algebra subspace 
SL(N, R)/SO(N) corresponding to a symmetric space of euclidean type (zero curvature) 
(cf. section ITT1|) and obtained by removing the set of real, antisymmetric and traceless 
matrices from the algebra SL(N, R). 

By performing a similar analysis for (3 =2, 4 we obtain the following general result: The 
gaussian ensembles labelled by /3—2, 4 consist of hermitean square matrices belonging to 
algebra subspaces SL(N, C)/SU(N) and SU*(2N)/USp(2N), respectively. From Ta- 
ble El we see that these three symmetric spaces correspond to algebra subspaces in Cartan 
classes A, AI and AIL The integration manifolds of the circular ensembles are the positive 
curvature coset spaces corresponding to the same Cartan classes. 

Also the chiral ensembles used in field theory correspond to algebra subspaces. They are 
identified respectively with SO(p, q)/(SO(p) <g> SO(q)), SU(p, q)/(SU(p) <g> SU(q)), or 
USp(p, q)/(USp(p) ® USp(q)) (in this case p, q have to be even). We will not discuss 
them further here (for details we refer to £Q). 

Example: The transfer matrix ensembles appear in the theory of quantum transport, in 
the random matrix theory description of so called quantum wires. In these pages we will 
only discuss the part of the theory which is relevant for our purpose, the study of the 
mapping between random matrix theory and symmetric spaces. 

The natural theoretical framework for describing mesoscopic systems is the Landauer the- 
ory [HS]- Within this approach Fisher and Lee |M] proposed the following expression for 
the conductance in a two-probe geometry (a finite disordered section of wire to which 
current is supplied by two semi-infinite ordered leads): 



H -> UHU" 1 + H' 



(19.8) 





(19.9) 



n 



56 



where t is the N x N transmission matrix of the conductor (see eq. (|19.11|) below) , N is the 
number of scattering channels at the Fermi level and T 1; T 2 • • - TV are the eigenvalues of 
the matrix tv . The Tj's are usually referred to as transmission eigenvalues. The constants 
e and h denote the electronic charge and Planck's constant, respectively. 

While the scattering matrix S relates the incoming wave amplitudes I, I' to the outgoing 
wave amplitudes O, O' (see eq. (119. 2j) ). the transfer matrix M relates the wave amplitudes 
in the left lead to those in the right lead: 



The transfer matrix formalism is more appropriate for description of Id systems than the 
scattering matrix formalism. This is due to the multiplicative property of the transfer 
matrix, as an infinitesimal slice is added to the quantum wire. 

Following a standard notation |2SJ EE] the scattering matrix has the following block struc- 
ture 




(19.11) 



where r, r' , t, t' are N x N reflection and transmission matrices. The unitarity of S implies 
that the four matrices tv , t't'\ 1 — rr\ 1 — rV' have the same set of eigenvalues 7\, T/v 
(0 < Ti < 1). The parameters 9 A, are related to the transmission eigenvalues by 



1 - T 

A l = (19-13) 



The Aj are non-negative. In terms of these, M can be parametrized as 



The Xi are the eigenvalues of the matrix 



Q = -(M f M + (M^M)- 1 - 2) (19.12) 
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where u, u', v, v' are unitary N x N matrices (related by complex conjugation: vl = u*, 
v' = v* if M G Sp(2N, R) or M G SO* (AN), see below) and A = diag(Ai, X N ). In case 
spin-rotation symmetry is broken, the number of degrees of freedom in (jl9.14j) is doubled 
and the matrix elements are real quaternions. 

Transfer matrices are strongly constrained by various physical requirements. Flux conser- 
vation, presence or absence of time-reversal symmetry, and presence or absence of spin- 
rotation symmetry lead to conditions on the transfer matrix. These conditions determine 
the group G to which M belongs. For example, flux conservation leads to the following 
condition on M jHEj 



i.e. M preserves the (2N x 2N) metric £ 2 . This means that M belongs to the pseudo- 
unitary group SU (N, N) (M has to be continuously connected to the unit matrix so we 
take the connected component of U(N,N)). For a detailed discussion of all the above 
constraints we refer the reader to pQ and references therein. 

Using the parametrization in eq. f)19. 14j) it is easy to check that rotating M by a matrix W G 
SU(N) x SU(N) x U(l) gives a new transfer matrix M' = WMW' 1 = U'TV with the same 
matrix V and therefore the same physical degrees of freedom {Ai, Xn}- This means that 
the matrix T belongs to a coset space G/K, where M G G and W G K. Following in each 
case a similar reasoning, one obtains three ensembles of transfer matrices corresponding to 
different physical symmetries: Sp(2N,R)/U(N), SU(N,N)/SU(N) x SU(N) x U(l) and 
SO*(4:N)/U(2N). These are symmetric spaces of negative curvature, as is also evident 
from Table 121 They correspond to Cartan classes CI, AIII, and Dill-even, respectively. 

In the three examples above, it was evident that the physical degrees of freedom in RMT 
correspond to points in the symmetric coset spaces in Cartan's classification. To every 
symmetric space a random matrix ensemble is associated, and many of these have been 
directly applied in some physical context (see pQ for more details). The ensembles listed 
in Table |U are examples of such RMT's. 

In section El we identified the random matrix eigenvalues with the radial coordinates on 
the symmetric space. Their probability distribution is determined by the Jacobian of the 
transformation to spherical coordinates. This Jacobian was given in eq. (jl8.17j) . It its de- 
termined by the curvature of the underlying symmetric space, and in addition by the roots 
and root multiplicities m a of the restricted root lattice (cf. eq. (|18.13jl ). As a consequence, 
the form and strength of the characteristic eigenvalue repulsion in RMT is determined by 
the root lattice of the underlying symmetric space! What's more, in the following examples 





(19.15) 
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we will also identify the Dyson and boundary indices of the random matrix ensembles as 
determined by these multiplicities, and we will see that the m a determine the classical 
orthogonal polynomials related to each matrix model. 

Example: The Jacobian associated to the chiral random matrix ensembles 

oc n w - \t n (w-ie) 

i<j k 

(where v is the number of zero modes) corresponds in the most general case to a root 
lattice with all types of roots. The chiral ensembles are algebra subspaces (this can be 
deduced also from the block-structure of chiral matrices). Like for the gaussian ensembles 
they are associated with the restricted root lattice of the curved symmetric spaces in the 
same triplet (since we have not defined any such root lattice for flat symmetric spaces). 
The restricted root lattices for the chiral ensembles are of type B q or D q for (3 = 1 and BC q 
or C q for j3 — 2, 4. The positive roots are as follows: for B q {e^, ± ej}, for D q {e^ ± ej}, 
for BC q {ej, ± ej, 2e;}, for C q {e,; ± ej, 2e{} (i ^ j always). Using the root multiplicities 
m = p, nil = P — 1, rn s = P\p — q\ = P v (cf. the values in Table |2J) we see that the 
Jacobian is of type ()18.23j) . It can be rewritten 

J {0) ({^})~n\^-^fl[\^\ a , P = m , a = m s + mi (19.17) 

i<j k 

From this Jacobian it is evident that in addition to the usual repulsion between differ- 
ent eigenvalues, Aj also repels its mirror image — Aj, and the eigenvalues are no longer 
translationally invariant. This kind of ensembles are therefore called boundary random 
matrix theories and a is called a boundary index. The boundary random matrix theories 
(BRMT) include chiral and normal-superconducting (NS) ensembles. We have given just 
one example, but one can make the same identifications in all the RMT's. 

Known symmetries of the random matrix ensembles can be understood in terms of the 
symmetries of the associated restricted root lattice. In particular, ensembles of A n type are 
characterized by translational invariance of the eigenvalues. This translational symmetry 
is seen to originate in the root lattice: all the restricted roots of the A n lattice are of the 
form (ej — ej). For all the other types of restricted root lattices (B n , C n , D n and BC n ) this 
invariance is broken and substituted by a new Z 2 symmetry giving rise to the reflection 
symmetry of the eigenvalues discussed above. 

As we saw in the introduction to random matrix ensembles, orthogonal polynomials are 
an important tool for the exact calculation of eigenvalue correlation functions in RMT. 
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Interesting in our framework is that the parameters which define the polynomials can be 
explicitly related to the multiplicities of short and long roots of the underlying symmetric 
space and thus, by the identification in Table EJ with the boundary universality indices of 
the BRMT. The relations are the following: 

Laguerre polynomials: 



L w (x) 



x x e x d n 
n! dx n 



'x n+x e 



X 



m. 



(x > 0) 
mi — I 



(19.18) 



Jacobi polynomials: 



(-l) n (1 - x)~ c d n ( {l + x) n+p 



v ' 2™ft! (1 + x)p dx n V(l-x) 



-n— a 



m s + mi — 1 mi — 1 
p = ^^ — , a = — (19.19) 



We see that A and p have the same expression in terms of m s and mi. Thus the BRMT's 
corresponding to Laguerre and Jacobi polynomials with the same A = p indices belong to 
the same triplet in the classification of Table El They are respectively the zero curvature 
(Laguerre) and positive curvature (Jacobi) elements of the triplet. This explains the fact 
that scaled correlators are the same for Laguerre and Jacobi ensembles of the same j3 near 
the boundary (so called "weak universality") pT7j . 

As a consequence of the identification of symmetric spaces and random matrix ensembles, 
the physical systems corresponding to random matrix ensembles can be organized into 
universality classes. In Table El random matrix ensembles with known physical applications 
are listed in the columns labelled X + , X° and X~ and correspond to symmetric spaces 
of positive, zero and negative curvature, respectively. Extending the notation used in the 
applications of chiral random matrices in QCD, where v is the winding number, we set 
v = p — q. The notation is C for circular, G for gaussian 10 , x for chiral, B for Bogoliubov-de 

10 Strictly speaking, these two groups are the scattering matrix and hamiltonian Wigner-Dyson ensem- 
bles, of which the latter have come to be referred to as "gaussian ensembles" due to the most common 
choice of random matrix potential, and the former arc called circular because the eigenvalues of the unitary 
matrices lie on the unit circle. 
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Table 3: Irreducible symmetric spaces and some of their random matrix theory realizations. 



Restricted 
root space 


Cartan 
class 


G/K (G) 


G /K (G^/G) 


m 


mi 


m s 


x + 


x^0 


X 


/1JV— 1 

A N -i 

A N -1 


A 

AI 
All 


SU{N) 
SO(N) 
SU(2N) 
USp(2N) 


SL(N,C) 
SU(N) 

SL{N,R) 
SO(N) 

SU*(2N) 

USp(2N) 


9 
1 

4 


n 





n 




c + 

^2,0,0 
^1,0,0 
W,0,0 


no 

^2,0,0 
pO 

^1,0,0 
pO 

^4,0,0 


2,0,0 
^1,0,0 
4,0,0 


IjG/, (v>a) 
w q \f ^ h t 


AIII 


£>U (p+q) 


° u (p,q) 


2 


l 


2u 




X2,l,2u 




C q (p=q) 


SU(p)xSU(q)xU(l) 


SU{p)xSU{q)xU{l) 


o+ 

2,1,0 


T~ 

2,1,0 


II AT 


B 


SO(2N + 1) 


SO(2N+l,C) 
SO(2N+l) 


2 


o 


2 




pO 

2,0,2 




Cn 
Cn 


C 
CI 


USp(2N) 

USp(2N) 


Sp(2N,C) 
USp(2N) 
Sp(2N,R) 


2 
1 


2 
1 






R+ 

D 2,2,0 

B ti,o 


R° 

r> 2,2,0 

R° 

1,1,0 


^2,2,0 
^1,1,0 


SU(N)xU(l) 


SU(N)xU(l) 


BC q (p>q) 


CII 


USp{2p+2q) 


USp(2p,2q) 


4 


3 






X4,3,4^ 




C q {p=q) 


USp{2p)xUSp{2q) 


USp{2p)xUSp{2q) 




T 4,3,0 


D N 
Cn 
BC N 


1 — 1 1 — 1 
1 — 1 1 — 1 
1 — 1 1 — 1 

Q Q Q 


SO{2N) 

SO(AN) 


SO(2N,C) 
SO(2N) 
SO* (AN) 


2 
4 
4 



1 
1 




4 


R+ 
r> 2,0,0 

R+ 

4,1,0 


r> 2,0,0 

R° 

4,1,0 

pO 

4,1,4 


^2,0,0 
^4,1,0 


SU(2N)xU(l) 
SO(AN+2) 


SU(2N)xU(l) 
SO* (AN +2) 


SU(2N+l)xU(l) 


SU(2N+l)xU(l) 






B q (p>q) 


BDI 


SO(p+q) 


SO(p,q) 


1 





V 








D q {p=q) 


SO(p)xSO(q) 


SO(p)xSO(q) 




T l,0,0 



Gennes, P for p-wave, T for transfer matrix and S for S-matrix ensembles 11 . The upper 
indices indicate the curvature, while the lower indices correspond to the multiplicities of 
the restricted roots characterizing the spaces with non-zero curvature. To the euclidean 
type spaces X° ~ G°/K, where the non-semisimple group G° is the semidirect product 
K Cg) P, we associate the root multiplicities of the algebra G = K © P. Note that the 
ensembles that have been given are the ones to which we have found explicit reference in 
the literature (see P). In principle all the empty boxes could be filled too. 

In Table 0] we list some identifications made between matrix ensembles and symmetric 
spaces. We have not discussed the Coulomb gas analogy, but most of the other entries we 
have or will mention. For a more thorough discussion and references we refer to pQ. 

11 For a more detailed discussion of all these ensembles and their relation to symmetric spaces we refer 
to pOj and references therein. 
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Table 4: The correspondence between random matrix ensembles and symmetric spaces 



Random Matrix Theories (RMT) 


Symmetric Spaces (SS) 


circular or scattering ensembles 


positive curvature spaces 


gaussian or hamiltonian ensembles 


zero curvature spaces 


transfer matrix ensembles 


negative curvature spaces 


random matrix eigenvalues 


radial coordinates 


probability distribution of eigenvalues 


Jacobian of transformation to radial coordinates 


Fokker-Planck equation 


radial Laplace-Beltrami equation 


Coulomb gas analogy 


Brownian motion on the symmetric space 


ensemble indices 


root multiplicities 


Dyson index (3 


multiplicity of ordinary roots {(3 = m ) 


boundary index a = f3(y + 1) — 1 


multiplicity of short and long roots (a = m s + m{) 


translationally invariant ensembles 


SS with root lattice of type A n 


boundary matrix ensembles 


SS with root lattices of type B n , C n , D n or BC n 


pair interaction between eigenvalues 


ordinary roots 



■* L *"—L— 

1 o 

Figure 8: Disordered wire of length L\ to which a segment of length L is added. This 
scaling operation leads to a Brownian motion of the transmission eigenvalues determining 
the conductance. Taken from ref. [TSJ. Used with permission. 

20 Solution of the DMPK equation 

One of the main problems in the RMT description of quantum wires is determining the 
probability distribution of the {Aj} variables appearing in the parametrization of the trans- 
fer matrix, eq. (|19.14|) . This gives access to the transmission eigenvalues and through the 
Landauer-Lee-Fisher formula, the main observable, the conductance G as a function of 
the length L of the quantum wire. To this end, Dorokhov |3S] and Mello, Pereyra and 
Kumar (35] derived a scaling equation which expresses the dependence of the probability 
distribution P({Aj}, L) of {X{\ as a function of L. This was done by considering the change 
in P({Xi},L) after addition of a thin slice Lq to the wire under certain assumptions (cf. 
FigEj). 
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The resulting equation reads 



dP 

ds 



DP 



(20.1) 



where s is the wire length measured in units of the mean free path /: s = L/l, and D can 
be written in terms of the {Aj} as follows: 



with 7 = (3N + 2-/3. /? is the symmetry index of the ensemble of scattering matrices, in 
analogy with the well-known Wigner-Dyson classification, and Jp(X) = J{{\}) is given 



The solution of this equation was not at all immediate. What's more, the appearance 
of Jp{{\i}) in the DMPK equation is due to the fact that the authors tried to mimick 
the Wigner-Dyson ensembles. However, this is quite misleading from the point of view of 
symmetric spaces, as the Wigner-Dyson ensembles have nothing to do with the transfor- 
mation from the transfer matrix ensembles to the space of the {A;}. As we will see, the 
final solution of this equation, which is exact for (3 = 2 and approximate for (3 = 1, 4, 
relies on the relationship of transfer matrices with symmetric spaces of negative curvature. 

The exact solution of the DMPK equation in the (3 = 2 case was first obtained in a 
remarkable paper [HHI by Beenakker and Rejaei. Here we review their derivation in a 
slightly different language, trying to stress the symmetric space origin of their result. 

The starting point is a mapping of the coordinates Aj onto the radial coordinates Xi on the 
symmetric space: 




(20.2) 



by 



(20.3) 



Xi = sinh 2 Xj 



(20.4) 



The DMPK equation can then be rewritten as 
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BP 1 

°— = -J{x)X B ,r\x)P (20.5) 



where 



J({%i}) = JJ|sinh 2 x i -sinh 2 x J -|^JJ|sinh2x fc | (20.6) 

i<j k 

is the Jacobian and A' B is the radial part of the Laplace-Beltrami operator on the underly- 
ing symmetric space! As a consequence, we can identify the differential equation with the 
equation for free diffusion on the symmetric space as a function of dimensionless "time" 
s = L/l where / is the mean free path related to diffusion in the quantum wire. 

At this point one can follow two roads. On the first one, which was the one followed 
by Beenakker and Rejaei, one maps the DMPK equation onto a Schrodinger equation in 
imaginary time. This requires the substitution 



P({x,}, S ) = JH{x t })q({x t },s) (20.7) 
A straightforward calculation shows that the DMPK equation then takes the form 



= (H-U)m (20.8) 

OS 



where U is a constant and 7i is a Hamiltonian of the form 



H = - - T (— + sinh- 2 (2x^ + ^ ~ 2) V Sinh2 ^ + Sinh2 ^ (20 9) 
n 27^\ar? + K ] ) 2 7 ^(cosh(2x i )-cosh(2a; J )) 2 1 U ' JJ 



At this point the main goal has already been reached: it is easy to see that if (3 = 2 the 
equation decouples: 



i o 2 i 

27 dx 2 27 smh 2 2x 
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The resulting equation was solved in jSH] using Greens functions. 



The second approach to a solution of eq. (|20.5j) also relies on the underlying symmetric 
space structure. Again, the starting point is the identification of the DMPK operator D 
and the radial part of the Laplace-Beltrami operator A' B on the underlying symmetric 
space: 



— = DP = —BP (20.11) 

ds 2 7 v ; 

where A' B = J- l BJ. 

As a consequence of this identification, if is an eigenfunction of A' B with eigenvalue 

k 2 , then J(x)§k{x) will be an eigenfunction of the DMPK operator with eigenvalue k 2 J (27). 
Then the properties of the eigenfunctions of the operator, which are zonal spherical 
functions, can be used to derive the exact solution of the DMPK equation for (3 = 2 
(|30j. see also |H]). This was done by Caselle. As the details of this calculation are quite 
technical, we will not review it here. The exact solution for the probability distribution is 
in both cases 



P({x n },s) = C(s) J^(sinh 2 Xj — sinh 2 xi) J^J(sinh2xj 



i<j 



x Det 



dk exp 



k 2 s 



tanh 



717 



.2m- 1 



AN I — v 2 . Pi (i ^i)(cosh2x„ 



(20.12) 



where P u {z) denotes the Legendre functions of the first kind. The second way of solving this 
equation constitutes a non-trivial consistency check on the solution obtained by Beenakker 
and Rejaei. 

The power of the description in terms of symmetric spaces becomes evident in the (3 — 1 
and P = 4 cases, in which the interaction between the eigenvalues does not vanish and the 
first approach discussed above does not apply. On the contrary, the description in terms 
of zonal spherical functions also holds in these two cases. Even though for (3 ^ 2 one does 
not know the explicit form of the zonal spherical functions, one can use an asymptotic 
expansion due to Harish-Chandra to get asymptotic solutions. The asymptotic results 
derived by Caselle can be found in 
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21 Relation to Calogero— Sutherland models 



An important role in our analysis is played by the class of integrable models known as 
Calogero-Sutherland (CS) models, which turn out to be deeply related to the theory of 
symmetric spaces. These models describe n particles in one dimension, identified by their 
coordinates q l , ...,q n and interacting (at least in the simplest version of the models) through 
a pair potential v(q l — qi). The Hamiltonian of such a system is given by 



where the coordinate q is q = (q 1 ,...,q n ), pi,...,p n are the particle momenta, and the 
particle mass is set to unity. In eq. (J21.1|) R + is the subsystem of positive roots of the root 
system R = {a 1 , a u } related to a specific simple Lie algebra or symmetric space, and n 
is the dimension of the maximal abelian subalgebra H . The components of the positive 
root a = a k G R + are a\, a*. The number of positive roots is u/2, where v is the total 
number of roots. In general, the coupling constants g a are the same for equivalent roots, 
namely those that are connected with each other by transformations of the Weyl group W 
of the root system. 

Several realizations of the potential v(q a ) have been studied in the literature (for detailed 
discussions see the review by Olshanetsky and Perelomov, ref. 0). We will be concerned 
with the following types of potentials: 



H 



i=l a£R+ 




(21.1) 



«i(0 

MO 

viii(0 



r 2 

sinh -2 £ 
sin" 2 £ 



(21.2) 



Example: The CS model corresponding to a C n root lattice is 
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^ = 1 y ^ 2 | y 9? \ ( 9 ° i 9 ° \ (213) 

2j-^d{q i ) 2 ^sinh 2 (2g i ) j-;. \smh 2 {q i - q j ) sinh 2 (g i + <f) J 

We remind the reader that the C n root system consists of root vectors {±2ej, dbe* ± ej, i ^ 
j}. The arguments of the sinh-function are q% = q-a k = (q 1 , q n )-(ei±ej) = q l ±qi (i < j) 
where a k is an ordinary positive root of the root lattice R, and (q l , q n ) ■ (±2ej) = ±2g* 
for the long roots. Note the different coupling constants for the ordinary and long roots. 

Under rather general conditions (see [H] for a detailed discussion) these models are com- 
pletely integrable, in the sense that they possess n commuting integrals of motion. As we 
will now see, the Calogero-Sutherland Hamiltonians can then be mapped onto the radial 
parts of the Laplace-Beltrami operators on the relevant symmetric spaces. The mapping 
is 



U = j\(q)-(M B ±p 2 ),r^(q) (+ for II, for III, p = for I) (21.4) 

where p is the vector defined in ()18.26|) and J(q) is the Jacobian for transformation to 
radial coordinates on the SS: 



rL eR +[<? c . 
= rL G R+[ sinh (<? c 

IW[sin(<f) 



)] ma II (21.5) 



Olshanetsky and Perelomov proved that this happens if and only if the coupling constants 
g a in H, take the following root values 



_ m a (m a + 2m 2a - 2)\a\ 



9l = ° V a ' g 71 ' (21-6) 

where m a is the multiplicity of the root a and |a| is its length. 

A number of results can be obtained for the corresponding quantum systems merely by 
using the theory of symmetric spaces. Due to equation (J21.4j) . once we know the zonal 
spherical functions, we can solve the Schrodinger equation. A detailed collection of results 
pertaining to spectra, wave functions, and integral representations of wave functions can 
be found in the original article jHj- 
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The equation f|2U.9|) in the solution of the DMPK equation discussed in the previous section 
can be recast in a slightly different form, thus completing the chain of identifications 

DMPK equation — symmetric space — Calogero-Sutherland model 
By using simple identities for hyperbolic functions, the Hamiltonian in (|20.9J) becomes |4*2*j 



n = V" ( 1 d 2 9? \ y> ( 9o 9 2 o \ (217) 

^\ 2dx 2 sinh 2 (2xi)/ ^ Vsinh 2 (x; - x, ) sinh 2 (xi + xA J 

where gf = —1/2, g 2 Q = (3{(3 — 2)/4 and c is an irrelevant constant. This Hamiltonian, 
apart from an overall factor I/7 and the constant c, exactly coincides with the Calogero- 
Sutherland Hamiltonian (|2Qj) corresponding to a root lattice R = {±2x,,, ±X{ ± Xj, i 7^ j} 
of type C n with root multiplicities m — [3, mi — 1. The values of the coupling constants 
g Q , gi are exactly the root values given in eq. (|21.fij) for which the transformation from TC 
into A' B is possible. 



22 Concluding remarks 



In these lectures we have given an elementary introduction to random matrix theory, as well 
as discussed some basic concepts in the theory of symmetric spaces. The scope of these 
discussions were to show, through numerous examples, that hermitean random matrix 
ensembles and the symmetric coset spaces based on simple Lie groups and classified by 
Cartan, are the same mathematical objects. By making this identification, new results can 
be obtained in random matrix theory by applying what is known in mathematics about 
such manifolds. We gave a few examples of this in the last lecture. 

In the first lecture we discussed a number of systems that are successfully described by 
random matrix theories. (Due to the limited expertise of the author, only some major 
applications in physics were discussed, however, there are also applications for instance 
in mathematics and biology.) We saw that random matrix theory gives a statistical de- 
scription of complex eigenvalue spectra of quantum operators in chaotic or many-body 
systems. By taking the double scaling limit in which the spectrum becomes dense and 
the eigenvalues are simultaneously rescaled on the scale of the average level spacing, one 
can separate out the universal part of the spectral behavior from the system-dependent 
one. We briefly discussed some typical spectral observables obtainable from the correlation 
functions. A large number of systems (some of which were discussed in the introduction) 
exhibit the generic spectral features arising in random matrix theory. 
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One such system is represented by quarks in a random gauge field background. We illus- 
trated this by showing a few graphs of computer simulations of Q CD-like gauge theories 
on a space-time lattice. In these the spectrum of the Dirac operator was studied and com- 
pared to analytical predictions from random matrix theory. Of course, there is a wealth of 
experimental and numerical results also on the spectra of nuclei, atoms, molecules, elas- 
tomechanical systems, microwave cavities, mesoscopic systems, etc., that we have not even 
mentioned and that could be given as examples of the type of spectral behavior that is 
predicted by random matrix theory. It is important to note that RMT behavior in spectra 
can also be partial, depending on circumstances not always fully understood (see [2] for a 
discussion). 

The typical eigenvalue repulsion in chaotic systems described by RMT can be traced to 
the geometrical properties of the corresponding symmetric space manifold. In particular 
we discussed how curvature, type of roots, and root multiplicities of the restricted root 
lattice on the symmetric manifold determine the exact form of these correlations. This 
follows from the general theory of coordinate systems on symmetric spaces and from the 
identification of the radial coordinates on the symmetric space with the set of random 
matrix eigenvalues. We made a number of further identifications between matrix ensembles 
and symmetric spaces and summarized the results in Table |U 

We also explained how symmetric spaces corresponding to compact symmetric subgroups 
are classified. All possible involutive automorphisms of the compact real form of an algebra 
leads, in a natural way, to a classification based on root lattices. This defines the univer- 
sality classes of the corresponding random matrix ensembles, and leads to a new scheme of 
classification dictated by the strict properties of root systems belonging to simple complex 
Lie algebras. 

The differential equation for quantum wires called the DMPK equation describes the evolu- 
tion with increasing wire length of the joint probability distribution of a set of parameters 
simply related to the transmission eigenvalues. In our discussion we concluded that it 
essentially corresponds to the equation for free diffusion on the symmetric space, and we 
discussed its solution in terms of zonal spherical functions for all three values of the Dyson 
index. The zonal spherical functions are known in the theory of symmetric spaces, where 
they play the role of eigenfunctions for the radial part of the Laplace-Beltrami operator 
on the symmetric manifold. 

Lastly, we discussed work by Olshanetsky and Perelomov on the integrability of a class of 
one-dimensional Calogero-Sutherland models. Such models become integrable at certain 
"root values" of the coupling constants, and the integrals of motion are then given by the 
Casimir operators related to the Lie algebra or symmetric space underlying the model. 
This leads to a set of exact results concerning spectra and wave functions of quantum 
systems. 
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It is certainly of interest to further investigate the possibility of applying known results 
on symmetric spaces to the corresponding ensembles of random matrices. An important 
research direction might be trying to fit also non-hermitean ensembles into a similar setting. 
This would bring us outside Cartan's scheme. In spite of the recent activity in the field of 
non-hermitean random matrix ensembles, the author is not aware of any research effort in 
this direction. There are also various types of extensions of Calogero-Sutherland models 
that could be explored in this spirit. 
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